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Summary—lIn the author's opinion, the 


“unit process ” 


‘ 


concept is completely against the 


logic of chemical engineering which pays no particular attention to the kind of atoms exchanged. 
Nevertheless, chemical changes cannot be left outside the field of that science. 
In an attempt towards a general classification of all unit operations, some can be put together 


in accordance with their mechanical character. 
momentum transfer, energy transfer, mass transfer, 


common character of exchange processes : 
intramolecular transfer. 


It is possible to classify all the others by the 


The proposed classification includes in the field of Chemical Engineering Science all the 
transformations carried out not only in the chemical industry but also in all the process 


industries. 


Résumé 


L’auteur se préoccupe de faire rentrer les transformations chimiques dans la systéma- 


tique du Gente Cuimigue, car il est peu satisfait du concept américain d’unit operation qui 
classe ces transformations suivant la nature des atomes échangés. 
Il est conduit & proposer une systématique générale englobant toutes les opérations fonda- 


mentales. 


Aprés avoir séparé un certain nombre d’opérations 4 caractére purement mécanique, 


l'auteur est conduit 4 classer toutes les autres sous le nom générique de phénoménes d’échanges : 

échanges d’impulsion, échanges d’énergie, échanges de matiére, échanges intramoléculaires. 
Cette classification permet de faire rentrer sous la discipline du Gente Curmigve toutes les 

opérations intéressant non seulement l'industrie chimique, mais; toutes les industries de trans- 


formation. 


We think that everybody will now agree, that 
Chemical Engineering is a true science : it is one 
of the many branches of the physical sciences, 
exactly like electronics, physical chemistry or 
Chemical engineering has its 
to conceive, to design, to 


hydrodynamics. 
own definite object : 
build and to operate the equipment with which 
any kind of chemical reaction is to be performed 
on the industrial scale. 

The idea of a true science of chemical engineer- 
ing was certainly in the air before the end of the 
nineteenth century, in some papers of G. E. 
Davis from Manchester, as well as in those of 
G. Lunce in Zirich or of H. Le Cuareier in 
Paris. But I think that the first clear exposition 





of the fundamental concept of chemical engineer- 
ing was made by Dr. A. D. Litt.e in his Report to 
the President of the Massachusetts Institute of 
Technology in 1915: 

““Any chemical process, on whatever scale 
conducted, may be resolved into a coordinate 
series of what may be termed ‘ Unit operations,’ 
as pulverizing, drying, roasting, crystallizing, 
filtering, evaporating, electrolyzing, and so 
on. The number of these unit operations is 
not large and relatively few of them are involved 
in any particular process. 

“The complexity of chemical engineering 
results from the variety of conditions as to 
temperature, pressure, etc. . . . under which 


* This was delivered as a lecture to the University of London at University College, London, in January 1953. 
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the unit operations must be carried out in 

different processes and from the limitations as 

to materials of construction and design of 
apparatus imposed by the physical and chemical 
character of the reacting substances.” 

There we have clearly defined the fundamental 
concept of our science: one basic operation 
carried out in different circumstances varying 
from the manufacture of one product to the 
manufacture of others. Before this concept was 
fully understood, the craftsman was trained by 
hard experience to know all the “rules of 
thumb” to deal with, say, a sulphuric acid 
sludge in a lead chamber plant but he would be 
left helpless if asked to deal with the sweet 
liquor in a sugar refinery. 

With the concept of unit operation, one has to 
understand that filtration is the kind of manipula- 
tion dealing with the mechanical separation 
through a porous medium of a fluid and the solid 
particles in suspension in it. Once the laws 


governing such a separation have been discovered 
and expressed mathematically in terms of general 
properties (density, viscosity, porosity, grist, 


etc. . . .), it suffices to give their special values 
in every particular case to all the numerical 
coefficients present in our mathematica! relations. 
The design of a filter, adapted to any particular 
problem, will be possible through simple calcula- 
tions. 

As Dr. Lirt.e said, the number of unit opera- 
tions is not very large. It has not increased much 
since 1915. But we have to face something quite 
strange : a scientific discipline intended to control 
all the past and future developments of the 
chemical industry rests completely upon a 
scientific concept which leaves out all the changes 
which are purely chemical by nature. Is it not 
inconvenient that the chemical changes, which 
will remain the essential foundation of the 
chemical industry, are to be left outside the field 
of chemical engineering ? 

A few of our American colleagues have tried 
to overcome that omission by the introduction of 
the additional concept of unit process. By unit 
process, they attempt to introduce the complete 
series of chemical changes like hydrogenation or 
reduction, oxydation, nitration, sulphonation, 





halogenation and so on. Personally, I found it 
very difficult to make our students pass from the 
very clear-cut concept of unit operation to the 
new concept of unit process. By the unescapable 
diversity of practical conditions under which they 
are to be performed in industry, the chemical 
changes defy any attempt at classification. What 
is there in common really between the nitration of 
toluene to make trinitrotoluene with a mixture 
of sulphuric and nitric acid and the vapour phase 
nitration of paraffins with nitrous fumes? Is 
there really something similar between the com- 
bustion of ammonia over a platinum catalyst and 
the cleavage of a double bond by aqueous per- 
manganate, both chemical changes belonging to 
the oxidation unit process ? 

In my opinion, it is useless to attempt to classify 
the chemical changes, by the nature of the atoms 
exchanged between different molecules; I con- 
sider definitely that the concept of unit operation 
is completely against the logic of chemical 
engineering. Our science pays no particular 
attention to the kind of atoms exchanged; it 
considers only the modalities of the chemical 
reactions. 

Nevertheless, it remains essential to introduce 
the various chemical changes in the methodology 
of our discipline : chemical engineering cannot 
leave chemical reactions outside its field of interest. 


Table 1. 


Mechanical unit operations 





Size alteration 
Reduction Enlargement 
Crushing-Grinding Pelleting— Briquetting 
Disintegration-Pulverisation Moulding-Sintering 
Cutting- Machining Extruding-Laiminating 
Pressing Flaking 
Heteregenous Systems 
Separation 
Sampling Screening 
Classification—Jigging 
Tabling 
Sedimentation-Decantation 
Elutriation- Filtration 
Centrifugation- Expression 
Flotation 
Separation (magnetic or 
electrostatic) 
Dust and mist collection 


Mixing 
Kneading 
Wetting Dispersing 
Spraying 
Atomisation 
Emulsification 
Homogenisation 
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In an attempt to overcome the difficulties of the 
unit process concept, I have tried to classify all 
the classical unit operations. 

I shall immediately set apart a class of unit 
operations which are of a purely mechanical 
character (Table 1): some are utilised when we 
need to perform an alteration of form or of size : 
crushing, grinding, pulverizing, briquetting, sin- 
tering, laminating, extruding, etc. ... We can 
also perform some other mechanical unit opera- 
tions when, dealing with heterogeneous systems, 
we want to separate or to mix various phases : 
filtration, sedimentation, centrifugation, mixing, 
blending, atomizing, etc. .. . 

For the remaining classical unit operations, like 
distillation, evaporation, cristallization, extrac- 
tion, and so on, I think that their nature presents 
identical characters and I propose to class all of 
them under the French name “ phenoménes 
d’échange ”’ which I shall try to translate by 
“exchange operations’ or by “ transfer opera- 
tions.”” It is necessary to note that in all these 
unit operations, our main interest lies in the 
kinetics of the transfer. Chemical engineering 
science is always obliged to investigate the 
migration of different molecular species through 
a surface of separation. 

The introduction of the new concept of “ trans- 
fer operations *’ appears purely as a matter of 
course when we stop considering chemical industry 
from the analytical side which made it possible to 
decompose any manufacture in series of unit 
operations. We must as well study the chemical 
industry in a synthetic way and try to find out 
what in it is always identical, no matter which 
products are considered. 

Any chemical manufacturing process is pri- 
marily a problem of traffic. For the production of 
any kind of chemical, one of the first problems 
lies in the introduction of the raw materials and 
the removal of the final products. We have also 
to introduce into our equipment some form of 
energy: usually heat is wanted or must be 
extracted from our plants. Flow of fluids, heat 
transfer are major chapters of chemical engineer- 
ing science. But it is obvious to me that here we 
have again some kind of physical “ transfer 
operation.” 





For these two kinds of physical transfer 
operations, we are led to utilise mathematical 
relations of identical form: the relations of 
similitude will introduce the dimensionless num- 
bers of Reyno.ps, of Nusseit, of PRaNDTL or 
Grasnor. To explain convection phenomena, 
we are led to see a kind of barrier, or “ transition 
layer’? and we introduce the concept of the 
“ film ” discontinuity at the surface of separation : 
it is the same film which can be rendered visible 
by some experiments on the flow of fluids. 

To interpret the diffusion of different kinds of 
molecules through the surface of separation be- 
tween two phases, we will again utilize the 
properties of this special “film.” It is perfectly 
logical to place together in the same chapter 
molecular transfer and momentum transfer or 
kinetic energy transfer: at the molecular scale, 
all these three exchange processes depend solely 
on individual contacts between real particles. 

Finally, I can class together under the common 
heading of physical transfer processes three 
different kinds of physical transfer : momentum, 
energy, molecules (Table 2). 


Table 2. 
Physical Exchange Processes 





Momentum Transfer 
Flow of fluids 

(gases and ideal liquids) 
Pumps—Fans—Compressors 
Non-Newtonian Systems 
Porous Beds—Suspensions 
Pneumatic transport 
Hydraulic transport 
Mechanical handling 


Energy Transfer 
Radiation 
Convection 
Heat exchangers 
Heat accumulators 
Thermophysics 
Heat engines 
Refrigerating engines 


Mass Transfer 

Single Systems 

Fusion—Solidification 
Sublimation 
Vaporization 
Liquefaction 


Complex Systems 
Distillation—Rectification 
Solvent Extraction 
Dissolution—Crystallization 
Gases—absorption and 

washing 
Solids absorbents 
Drying 
Conditioning of gases 





But any chemical reaction is also an exchange 
process by which we interchange different atoms 
between different molecules. But chemical 
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Table 3. 
Chemical Exchange Processes 





Fluid Systems 
— 


Reactors 


= 


Batch operation — 


Continuous reactor ———_ -- 


Tubular 


Multistage 


Countercurrent - 
Operation 
Cocurrent 


Isothermal 


Pseudo isothermal 


Systems with Solids 


Adiabatic —— - 4 


Converters 


4 
Batch operation 


Continuous converter ’ 
| 
Single Multistage Tubular 


Immobile 

. , Bed 
Semi-mobile J 
Countercurrent 


* | Operation 
Cocurrent | 





Fluidised Bed 


r- — — Converter 


Single 


Multistage 


Tubular 





engineering will pay no particular attention to 
the kind of atoms involved and remains interested 
only in the modalities of the transfer. As for all 
other transfer processes, chemical engineering will 
consider in the chemical changes, the driving 
force and the rate of exchange. All the chemical 
exchange processes will be classified according to 
the type of kinetics imposed by the circumstances 
under which they occur. These circumstances 
will vary from one reaction to the next but the 
nature of rate exchange will dictate the type of 
equipment required. Our classification of the 
chemical transfer processes will depend mainly 
on the different types of equipment: according 
to type of energy or mass transfer we are com- 
pelled to utilize different kinds of equipment. 
For some particular chemical reactions, we 
want to prevent any heat exchange with the 
outside ; we need perfect thermal insulation : we 
will attempt an adiabatic operation. On the 
contrary, in some other cases, we need to main- 
tain the temperature constant and we can achieve 
our isothermal transformation if we arrange 
to provide adequate heat transfer through the 
walls. Between these two extreme conditions of 





operation, practical necessities will introduce 
intermediate imperfect modes of transformation : 
if the wall alone can be kept at constant tem- 
perature, we operate a pseudo-isothermal trans- 
formation. 

About the mass transfer, the circumstances van 
also be extremely variable. In batch processes, we 
introduce in our equipment all the raw materials 
in the required sequence and after a certain 
delay we take away the final products. Chemical 
industry now utilises more and more continuous 
In the proper equip- 
ment, we feed at a steady rate the raw materials 
and obviously the character of our equipment 
will be dictated by the flow characteristics of our 
raw materials ; if we deal exclusively with fluids, 
the limitations of the equipment are less tight 
than if we have to deal with solids. 

I consider that these differences in the flow 
characteristics of the raw materials, are so 
important for chemical engineering that I 
will introduce a more precise terminology. I shall 
give the generic name of “reactor” to the 
equipment dealing exclusively with fluids either 
homogeneous or heterogeneous. I shall reserve 


processes of manufacture. 
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the name of “converter” to the equipment 
dealing with mixtures of solids and fluids. Ob- 
viously converters and reactors can be run either 
batchwise or continuously. They can also be run 
either adiabatically or isothermally, or pseudo- 
isothermally. 

For the continuous reactors, where we deal 
either with homogeneous or heterogeneous raw 
materials, we can consider first a type where they 
will stay long enough to give time for the chemical 
change to be achieved : 
single reactor. With a perfect stirrer, at every 
instant, we will find no variation of composition 
between different locations in our reactor. 

In many cases, we may be obliged to fractionate 
our reactor in different compartments, each of 
them performing a definite step in the complete 
chemical transformation : we then deal with a 
compound reactor or multistage reactor. This 
equipment works in a completely different way to 
the single reactor: it makes easier the isothermal 
operation and it will be utilized if the chemical 
reaction requires different temperatures for the 
different steps. If the number of compartments 
or stages grows without limit, we arrive finally 
at the tubular reactor ; in this type of equipment 
the fluids will flow at a constant rate and in every 
section we will reach a steady state of concentra- 
tions, velocities, temperatures, and so on. 

With the multistage or tubular reactors, we 
can vary at will the local concentration of our 
reactants, as we can choose various feed points 
for our different raw materials. These feed points 
can be located at the opposite ends of our equip- 
ment and we will work counter-current wise, as 
opposed to cocurrent operation when all our 
initial reactants are fed at the same end and the 
final products withdrawn at the other. 

If we deal with solid materials, we shall use 
converters which can be either adiabatic, iso- 
thermal or pseudo-isothermal. If our solid 
materials remain immobile, we have a fixed-bed 
converter : we can find many examples of fixed- 
bed converters in catalytic conversions. We can 
deal with single converters or with multistage 
converters in sulphur dioxide oxidation. In 
ammonia synthesis, we deal usually with tubular 
converters. 


we have a continuous 





But in many other cases, the solid phase will 
play its part in the chemical change, but usually 
its movement is very slow compared with the 
flow of fluids. Our bed will be a semi-immobile 
bed : gas producers, blast-furnaces, rotary kilns 
or dryers are typical examples of our semi-immo- 
bile bed converters. 

In scientific matters, some kind of classification 
is absolutely necessary to clarify our problems, 
but we must always remember that nature does 
not care for the limitations of our brains. Classi- 
fication, always, must remain flexible. With the 
fluidisation we have the best example of the 
transition step between the two 
extremes: reactors and converters. A fluidised 
bed equipment deals with solid particles and 
belongs legitimately to the converter class. But 
fluidized beds partake of the general properties 
of fluids and I cannot raise objections against 
fluidized bed converters being called also fluidized 
bed reactors. 


necessary 


It is not very long since chemical engineering 
has begun to investigate the characteristics of the 
different chemical exchange processes and to make 
the necessary distinctions between the different 
types of equipment best suited to different 
chemical reactions. 


It gives me a real pleasure to pay a well 
deserved tribute to Dr. Densicu who was the 
first, I think, to give the lead in this new develop- 


ment of chemical engineering. In a paper pub- 
lished in the Transactions of the Faraday Society 
in 1944,* he put forward the theory of the 
continuous multistage reactor, a theory for the 
nitration process based on researches carried out 
during the war in the Royal Ordnance factories. 
DenBiGu extended very soon his theory to poly- 
merization processes: the same type of theory 
deals adequately with changes involving very 
different kinds of chemical transformations. 

In the United States, chemical engineers are 
now devoting their attention to the theory of 
converters, and we begin to learn how to design 
with some accuracy either adiabatic or isothermal 
or pseudo-isothermal converters. 


- * Dension, K. G., Trans. Faraday Soc., 1944 40 


352-373 ; Faraday Soc. Discus3., 1947, p. 263. 
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I really hope that my personal attempt at a 
close comparison between momentum or kinetic 
exchanges and molecular or atomic exchanges 
will not remain purely formal. The concept of a 
converter or a reactor especially adapted to a 
particular class of exchange is not really limited 
to purely chemical exchanges. 

Let us look at the most usual equipment for 
mass exchanges between liquids and gases: the 
packed tower is usually employed for the physical 
absorption of soluble components, like ammonia 
or sulphur dioxyde in water. But it is also 
utilised when the absorption follows some chemi- 
cal reaction as in the case of carbon dioxide and 
caustic solutions. What could tell the differences 
between a packed tower and a continuous tubular 
reactor ? 

The bubble-cap plate tower is not basically 
different from Densicn’s multistage reactor, 


except as to the efficiency of the mixing between 





vapour and liquid. Perhaps this will open a new 
way of attack to a very old problem, still unre- 
solved. In distillation, we know with a gratifying 
accuracy the equilibrium relationships between 
liquid and vapour compositions but we know 
very little about the rate of approach towards 
the equilibrium in our rectifying columns, in spite 
of what we can hear or read about theoretical 
plates, height equivalent to theoretical plates, 
height of transfer unit, Murphree efficiencies or 
over all efficiencies. 

Perhaps, once the real identity between plate 
towers and multistage reactors has been accepted, 
it will be easier to find some way out of our 
present difficulties on the kinetics of distillation 
problems. Perhaps the applications of Densicn’s 
line of attack about by-passing losses might give 
a quantitative account of the departure from 
equilibrium between liquid and vapour when we 
mix them on the plates of our stills. 
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On a problem of diffusion and simultaneous gas flow 


W. Jost 
Institut fiir Physikalische Chemie Universitat Gottingen 


(Received 1 September 1953) 


Summary—A diffusion problem of spherical symmetry, involving gas flow, has been treated 


quantitatively. 


Résumé—L ‘auteur propose une solution quantitative pour un probléme particulier de diffusion 
& symétrie sphérique, accompagnée d’un écoulement gazeux. 


The following calculations were suggested by a 
paper presented by W. Morawtetz at the Duis- 
burg meeting of the Deutsche Bunsen-Gesell- 
schaft.* If the reaction 


8 FeS + 2S0, + Fe,O, + 2.55, (1) 


is carried out with a sphere of FeS in a stream of 
inert gas, containing SO,, the rate of reaction 
is determined by the diffusion through an outside 
porous layer of Fe,O,. The diffusion proceeds in 
an outward gas flow because for every 450, 
molecules which enter the reaction, 5 molecules 
of S, leave the interface FeS/Fe,O,, where the 
heterogeneous equilibrium according to (1) is 
assumed to be established. This diffusion problem 
can be treated without simplifications, as follows. 
Instead of (1) we use the general equation : 


S, +¥G, +> S, +4 G,. (1) 


Where a solid S, reacts with a gas G,, and yields 
the solid S, and another gas G,, the number of 
moles » of the newly formed gas may be greater 
or smaller than that of the consumed gas y, cor- 
responding to an outward or inward flow of gas. 

In Fig. 1 p denotes the radius, r is the variable 
radius of the inner sphere of solid S,, and r, is 
the constant radius of the outer porous layer of 
solid S,. The speed of gas flow in the direction 
of + p is v, depending on p, according to 


r 
U =U, * (2) 
p 


=v r,? 


a 
UP oo? 


if v = vy, for p=r,. The concentrations in the 
gas phase (within the porous material) are c, and 
¢, for G, and G,, c, + ¢, = ¢, = const. When 


* To appear in Z. Elektrochem. 





D is the effective diffusion coefficient for the 

porous material the diffusion flow is 

J, = D oe). 
dp 


dC, 


3. _ J): : 
J; “ 


(3) 


YY 


Fig. 1. Experimental set-up, used by Morawietz, 


schematically, and notations. 


The flow due to convection is 


and the total flow 
Ji 


J, 


I, + J,= (4 + €) %-% 
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The time-dependent diffusion equation for 
component 1 is 


2 2 ’ 2 > » 
1 _ _ dif J,) =D] 214.7%) _Yofo M9 6) 
dt de® prp] = p*® %%~ 


which may be put equal to zero for a quasi- 
stationary state. 
In this case c, depends on p only and we may 
write instead of (6) 
pic, + ¢, (2p —a) = 0 (7) 


where dashes denote differentiation with respect 
to p and we have introduced the length « 


a =v, 7,°/D (8) 


Equation (7) is to be solved for the following 
boundary conditions : 


forp=r:pJ,=—yv Jo, €1 Ce (9a) 
according to equilibrium (1), ‘ 

for p = 1, : Cg » 0 and (9b) 
(Jaro = Ia), = (9c) 


assuming that component 2 is carried away 
sufficiently quickly by the outside gas flow. From 
(9a) follows 


»|-0(% t) + ety 4 =-»| o(%) + este | 


and 


> 72 
p(*:) ~ eye : s + vegt 6%] /e—) (10) 


and from (9c), considering that c, = 0 for p=r, : 


> 2 : 
b (%) - + Cy0% b(*) (11) 
oP /rto P/ ro 


Putting c’ in eq. (7) equal to é we have 


dl —% 
me 2 2 ®, (12) 
dp p* 
and integrating between p = r, and p, we have 
ug tat ..? —2log?, (13) 
fo % Pp To 


£/£5 = (ro/p)® exp (a/r9) exp (— a/p). (18a) 





By a further integration 


? 


2 
¢, — ¢,° = | fap tafe 


. x [exp (« r,)exp(—a/p)—1] (14) 


Eq. (18) gives (2) (= €forp =r) expressed 
dp 


r 
by means of (“*) * Substituting this in eqs. 
dp 
(10) and (11) we have two equations for v, and 
(“*) , and with the values thus obtained we 
dp /ro 
may calculate c, from (14). With p = r and with 
« from eq. (8) eq. (18a) is written as follows 


(21) = (1) "exp (978 


and from eqs. (10) and (11) 


exp — rare) (15) 


d 
b(=*) “1 = pt Ug (Cy" + €y")/(m — v) (16) 


where c,’ and ¢,’ denote the concentrations for 
=r. Further, from eqs. (11) and (15) 


_ pf 21) % ’0) _ Me in 7 
r9=D| =) ft —exp( "0" | exp( -D)I (17) 


and with (16) 
HU (Cy" + Ca!) 


bv 


T =—- 
C2 % = 


Up To® 
it)] (18) 


Solving (18) for the exponentials and taking 
logarithms we finally have 


. _ oS 
Vo YT" ao 


“DD ey” +") v9) 


and 
) r r 

-/ ee, F (20) 
Tr BY +e") T—T 
Since ¢,"” = ¢y, the original concentration of 1 in 
the gas phase, and since c,” is to be assumed as 
given by the concentration according to the 
heterogeneous equilibrium at p = r for S, + v G, 
= S, +4 G,, ( =e," —¢,” is known, too. On 
the other hand, we have from (14) 
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C=¢," —¢,” 


are | 


where vu, is given by eq. (20). 


ss] 


{ and v, being 


known, (21) is an equation for (**) ° 
dp / ro 
We now calculate the outward flow of com- 


ponent 2 through the outer surface at p = ry. 
This is given, per unit area, by eq. (5b): 


) (c, being zero at p= T,.) (22) 
dp/ ro 


Substituting form eqs. (21) and (20) we find 


r r 
{ Driog Pe 794 
w(ey" + C2") | 


os Ug? 
role —ro)| exp| 


(22) 


From (28) with (19) 


Be” +e,’ 
_oDr wer” +") 
ro (r — To) (v — pm) Cy" /p (Cy" + C4") 


J, = (24) 
Eq. (24) assumes a particularly simple form for 
the special case (v — p)¢," < w(e," + ¢,"); in 
this case we write 


| pe,” + ve,’ 
# (ey" + €9") 


— Jog Liles” +62") +4(¥ —H)e4") +41» —H)eQ" 
[wley’ +9") +4(¥ —p)e2")] —4(v —p Der” 


1+(y—p)e4”, [2m(e" +¢,") +(v —p)e,"] 


1 —(v—p)eg” /[2p( ey” +¢9") +(» —)r2"] 


=... . 
2ye (€," + €2") + (v — p) ,” 


Thus finally 


{Dr 2u (¢,’ + ¢,’) 


te lr — re) Bale +e) + — Wee 


(27) 





The total flux through the outside surface of area 
47,7 is 


I= 4rr,? J,. (28) 
If d, is the density of phase 1, and d, density of 
phase 2, a change dr during time dt corresponds 
to a change dw in weight 


dw = 4m r* dr (d, — d,). (29) 

If in eq. (27) {, the concentration difference, is 
expressed in weight per cc., then (28) gives the 
total mass flow of component 2 out of the sphere. 
The rate of mass loss is obtained by multiplying 
this expression with («J/, — vM,)/uM,. Here 
M, and M, are the molecular weights of com- 
ponents 1 and 2. and use has been made of eq. 
(1). In the example quoted on the first page this 
expression is positive, but in other cases it might 
be negative as well. Thus we obtain from (28) 
and (29) : 


(30b) 


There are several possibilities for a practical 
evaluation of experiments. One may integrate 
(30b) between the limits t = 0(r = ry, w = wo) 
and t=7r (w=w,, r=r,), to give 


— = (re — 3) (dy —d,). 


(31) 
When w = w(t) has been measured as function 
of t, eq. (31) gives r = r(t), as function of t. In- 
serting these values of r in eq. (30a), where J, 
depends on r according to (27), we obtain an 
expression for dw/dt, still depending on an 
arbitrary parameter D. If D can be adjusted in 
such a way as to fit the experimental values over 
a sufficient range of t, the problem has been 
solved, provided that the value of D is reason- 
able, i.e., corresponding to gaseous diffusion in 
a porous medium. Then the underlying assump- 
tion, that diffusion is the rate determining step, 
may be assumed as verified. 
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If, on the other hand, we want to calculate the 
time necessary for a certain amount of reaction, 
we write eqs. (30a) and (30b) 


4m 1? (d, —4,) 7 = F(r) 


—vyit 
vM, 47 r? J, 


(32) 


F (r) = es 


uM, 
dt = 4m 7° (d, —4,) 
F (r) 
Integrating (32a) between t= O(r =r.) and 
t= 7(r =r.) we finally obtain for the reaction 
time +r 


rts {* r? (d, — d,°) = 
P(r) 


(32a) 


(33) 
ro 

The integration is to be performed graphically. 

The value of the effective diffusion coefficient D 

must either be taken from former experiments or 

be estimated appropriately. 

For a numerical evaluation of actual experi- 
ments we refer to Morawtetz’ original paper. His 
method differs slightly from that presented here. 
The above calculations present a justification of 
his attempt to take into account the gas flow by 
introducing independent diffusion coefficients for 
both components. 





NOTATION 


a = abbreviation for vgrs”/D, of dimension of a 
length. 
¢). Cg = concentrations of components 1 and 2 in the gas 
c’ = d,/dp = €. 
D = diffusion coeflicients [em? sec™"), 
G,, Gy = gaseous reactants, of concentrations c, and ¢y. 
1 = total flux of component 2 through outer sphere 
of area 4x 19°. 
J = flow vector (mass passing in positive r-direction 
per unit time perpendicularly through unit area) 
indices 1 and 2 referring to components 1 and 2. 


Ja - Gaia tee upper indices 1, 2 referring to 


‘ components 1, 2, respectively. 
J. = convection flow I a ’ 


M,, M, = molecular weights of components 1 and 2, res- 
pectively. 
,» = coefficients in reaction equation. 
p = radius, variable. 
r = variable radius of inner sphere. 
To = constant radius of outer sphere. 
S,, Sg = solid reactants. 
v = flow velocity. 
&=c’ = de,/dp. 


v 


The author would like to use this opportunity to correct 
a mistake overlooked in his book (Diffusion in Solids 
Liquids and Gases, Academic Press, New York, 1952) Eq. 
[1.237]on p. 54 should read 


v% To de 


D [d%e/dr® + (1/r) de/dr] — =0 


r ow 
one solution of which is ¢ = 8 r°°'°!?, with 8 arbitrary 
constant. 
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Hydroextraction VII: Residual moisture in whizzed cakes 
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Summary—The consolidation of fine particulate beds due to compressive stresses makes it unlikely 
that a general relationship can be found for the de-watering of such hydroextractor cakes by whizzing, 
though empirical linear functions can be used for the effect of rotational speed on residual moisture 


content. 
Résumé 


Il ne parait pas vraisemblable que l'on puisse établir une relation générale pour l’égouttage 


des gateaux par le turbinage, a cause du tassement des lits de particules fines prov qué par les efforts de 


compression. 


Cependant des équations linéaires empiriques peuvent donner l'effet de la vitesse de 


rotations par ’humidité residuaire :u cours de (urbinage. 


An important part of a hydroextraction cycle is 
the de-watering of the cake by “ whizzing” 
without liquid being fed to the basket. The 
following studies were to assess the possibility 
of using simple correlation functions of appropriate 
variables such as surface tension and angular 
velocity to describe the variation of residual 
moisture content in cakes of fine particulate 
material. The initial compressibility of the cakes 
was also studied by observing the dependence of 
the moisture content of submerged cakes on 
hydroextractor speed on first running at success- 
ively higher speeds. The latter is of interest in 
assessing the packing rearrangements occurring 
in submerged cakes. The highest compression 
attained under these conditions is also of signi- 
ficance in the relationship between moisture 
content and whizzing speed in cakes running 
dry. 


EXPERIMENTAL TECHNIQUE 


In procedure A the granular material (chalk, 
maize starch or kieselguhr) was made to a thin 
slurry and fed to a 9 in. diameter hydroextractor 
basket with the precautions necessary to induce 
the formation of a smooth cake of vertical 
surface [1]. The centrifuge was spun at 900-1000 
r.p.m., and thereafter for 15 minutes with the 
cake submerged under water. Without allowing 
the cake to run dry, the basket was stopped and 
samples of the cake taken out with an annular 
cutter of 1 in. diameter. Samples were cut from 
different circumferential positions and heights in 
the basket. Each sample was sliced into three 





approximately equal parts to represent the inner, 
centre and outer radial positions of the cake. The 
weighed samples were dried at 105°C to constant 
weight. 

Immediately after taking samples from the 
cake, it was run again submerged at a higher 
speed of rotation with the liquid inner surface 
being maintained at the same radius. After 5 
minutes running at the new speed, samples were 
taken and treated in the same way as above. 
The procedure was repeated over the speed range 
of the hydroextractor. 

In procedure B, after forming a cake in a 
hydroextractor basket, the speed was increased 
gradually up to the required value and maintained 
for 15 minutes with the cake still submerged. 
The liquid feeding was stopped and the cake was 
spun at a certain speed. Samples of complete 
cake thickness from mid-height were taken every 
minute, and their water content measured gravi- 
metrically. It was found that 12 minutes was a 
reliable time of spinning for reaching a constant 
moisture content at different speeds and surface 
tensions. 

In procedure C, after forming a cake and 
running it submerged at a certain speed, it was 
run dry (without liquid feed) for 12 minutes at 
the lowest speed, and samples were taken. The 
cake was then run at a higher speed, and sample 
taken, the procedure being repeated up to the 
highest speed (2800 r.p.m.). 

In procedure D, a series of small cakes was 
formed by suction in a 1} in. dia. cell [2] which 
was then mounted on the hydroextractor and 
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tested in the same way as mentioned in pro- 
cedures B, C. Air was sucked through other 
samples in the small cell and the residual moisture 
contents obtained for comparison with the whiz- 
zing data from procedure B. 
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Fig. 1. Moisture content distribution in hydroextractor 
cakes ; (a) along the depth of the basket : T-top; M— 
middle ; B—bottom ; (b) in a radial direction ; F-—front, 
C—centre ; B—back. 











Where possible the series were performed with 
maize starch, precipitated chalk and kieselguhr. 
To study the effect of surface tension Teepol 
was added to the formation and wash waters. 





EXPERIMENTAL Resucts: ror Wer Cakes 
The “ moisture content,” R, of a sample is 
defined and used hereafter as the ratio of the 
liquid to dry solid in the sample, expressed as a 
percentage. This value is directly proportional 
to the void ratio in a sample which is completely 
submerged. The voids available to flow in a 
tightly packed fine solid may differ appreciably 
from the value appropriate to the above ratio. 

Samples from either wet or spun-dry cakes 
gave reproducible moisture contents when taken 
at the same depth in the basket for any circum- 
ferential position. 

The moisture content for samples from pro- 
cedure A varied with the depth of basket, being 
always greater at the top. This variation was 
insignificant in case of starch cakes, whereas 
chalk and kieselguhr showed 10-15% higher 
moisture content at the top than at the bottom 
of the cake. The rate of change of moisture 
content with depth was higher in the top half of 
the basket, as shown in Fig. 1. 

The radial distribution of moisture content 
passed through maxima in the cases of chalk 
and starch with the value at the inner radius 
slightly less than that at the outer radius (Fig. 1) 
and the maximum about 5% higher than the 
value at the inner radius. Kieselguhr showed a 
gradual decrease of water content from inner to 
outer radius (Fig. 1). 

Chalk cakes formed under hydrostatic heads 
in the small cell did not show the type of variation 
noticed in cakes formed in the basket. The water 
content decreased 20°%, from front to back of the 
cake (Fig. 2, data: A). Samples formed by 
feeding to the spinning small cell and spun sub- 
merged gave similar moisture distributions such 
as the data B in Fig. 2. Agreement between the 
hydrostatic and spinning data was obtained from 
both chalk and starch cakes, showing that forma- 
tion method has no influence on the packing. The 
maxima obtained in the distributions for hydro- 
extractor cakes (Fig. 1) are thus indicative of 
packing characteristics peculiar to the hydro- 
extractor cake. 

The moisture content of samples from procedure 
A was independent of the speed at which a wet 
starch or kieselguhr cake was formed and washed. 
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The change was within 2°, when the centrifugal 
stresses were increased six times (Fig. 17). 
Chalk showed a gradual decrease of moisture 
content with spinning speed (Figs. 3 and 4). The 
water content after running at n’ = 2500 r.p.m. 














CONTENT 
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PERCENTAGE MOISTURE 








___ Front | 


= 
Fig. 2. Radial distribution of moisture content in cakes 
formed in the small cell; A: by feeding slurry to the 
by feeding slurry under hydrostatic 
suction head. 


spinning cell. B: 


was 20°, lower than after running at n’ 1000 
r.p.m. The rate of change was of the same order 
at the inner, centre and outer radial positions of 
the cake, with the same distribution mentioned 
above (Fig 1) holding at most of the speeds. In 
some cases, however, the moisture content at the 
outer radius of the cake was higher than the 
middle, as at speeds the 800 and 1600 r.p.m. 
(Fig. 3). The test on the cake represented by 
the dotted curve B, Fig. 4 was repeated after 
attaining the highest speed. The difference bb, 
represents the change due to altering the level of 
samplings, and the curve B, is the result of the 
repeated test. The condition attained after the 
original test was completed was hardly altered by 
repeating it, showing that the compression effect 
is not elastic. When Teepol was used to halve 
the surface tension of the formation liquid, the 
cake retained less water but the shape of curve 
remained the same (Fig. 4). 


Although smooth curves drawn 


may be 





through the experimental points on Fig. 3, 4, yet 
the abrupt changes cannot be ignored near the 
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Fig. 3. Variation of moisture content with speed for a 
submerged chalk cake ; samples from mid-height of the 
basket. 
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SPEED Wx IO" p.m. 
Fig. 4. Variation of moisture content with speed for 
submerged chalk cakes ; “through”’ samples from mid- 
height of the basket ; A, B, C : different cake thicknesses. 
B, : result of repeating test on cake B ; samples from top 
of the basket, 6, b; represent the difference due to change 
of sampling height. D : cake containing the same weight 
of chalk as B and formed from water whose surface tension 

was lowered by an addition of Tecpol. 
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speeds 1300 and 1800 r.p.m. On a semi-log plot 
(Fig. 5) these changes are clear. 


The relation 
R= R, — C. log, n’ (1) 


could represent the variation at different stages, 
with R, and C changing from stage to stage. It 
is most unlikely that equation (1) can be taken 
to represent a mean straight line through all the 
experimental points. The cake apparently suffers 
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Fig. 5. Semi-logarithmic plot of the results in Fig. 4. 


abrupt porosity changes caused by redistribution 
of the packed particles. 
The water saturating a chalk bed will exist as : 
(a) the equilibrium water content when sub- 
merged, i.e., water absorbed by the particles 
themselves, 
(b) water filling the interstices between the 
particles and removable by flow, 
(c) water held to the solid surface and not re- 
movable by flow. 
The void ratio significant in flow phenomena is 
only dependent on (b) and a value representing 
the items (a) and (c) should be subtracted from 
the total water content in measuring the void 
ratio by moisture content. However, if (a) and 
(c) are considered constant, which is a reasonable 





assumption, the curves in Fig. 7 would not be 
changed whether drawn with R or the quantity 
representing (b) as ordinate. 

The values (a), (c) are difficult to determine. If 
there is no evaporation during the process of 
spinning a cake dry, then the equilibrium value 
at the highest speed, R,, will be mainly (a) and 
(c). Fig. 6 is a plot of log (R — R,) against 
log nn’, and shows that the relation between 
void ratio and stress is not a simple power func- 
tion. There is an apparent anomaly in that the 
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Fig. 6. Logarithmic plot of the results in Fig. 4. 


curve for the cake formed from lower surface 
tension liquid is at higher ordinate than the 
others, leading to the conclusion that lowering 
the surface tension increases the porosity, which 
is contrary to the permeability result reported 
before. It is likely that during the dry running 
of the cake it is not only (b) that is affected, but 
(c), and perhaps (a) as well. 


EXPERIMENTAL Resutts ror Dry CAKEs 


Whether cakes are dried by sucking air through 
them or by spinning, a certain time has to elapse 
before the moisture content reaches a constant 
value for a specific extractive force. Although 
the first few moments are the most effective in 
the process, yet the path of the curve cannot be 
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due to packing changes 


Twelve minutes is required for equili- 


predicted 
(Fig. 7). 
brium to be reached in the cases represented in 
Fig. 7. 

When the cakes were run dry, the distribution 
of water content from top to bottom of the cake 
remained similar to that when wet, but the radial 


apparent 


distribution changes as shown in Fig. 8. The inner 
part was usually the driest in the cakes of kiesel- 
guhr and starch. Chalk reached constant dryness 
all through. With starch and kieselguhr a dry 
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Fig. 7. Variation in residual moisture content with time : 
A: chalk cake formed in hydroextractor, washed with 
water and spun at 800 r.p.m. B : 5 g. chalk cakes formed 
by suction in small cell, washed with water and spun 
dry at 800 r.p.m, C : chalk cake formed from water-Teepol 
solution, in hydroextractor and spun at 1000 ¢.p.m. D :5g 
chalk cakes formed by suction in a small cell, washed with 
water and dried by sucking air through the cake. E : Starch 
cake formed in hydroextractor, washed with water and 
spun at 900 r.p.m. F : 5 g chalk cakes formed from water- 
Teepol solution, by suction in the small cell and dried by 
sucking air through the cake. 





crust was observable at the inner radius of the 
cake. The material susceptible to increases in 
stress when running submerged, gave constant 
dryness across the cake, as with chalk, whereas 





the front layer forms a loose crust with starch and 
kieselguhr, which were materials unaffected by 
increased stresses when wet. 

Figs. 9, 10, 11, 12 and 13 represent tests on 
one chalk cake formed from 1000 g dry chalk 
which was made to a thin slurry and fed to a 
hydroextractor rotating at 1000 r.p.m., then spun 
wet at 2000 r.p.m. Fig. 1 shows the wet condition 
of this cake. The data in Figs. 9-13 were obtained 
by procedure C, Increasing the speed of rotation 
when the cake is dry causes a movement of water 
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Fig. 8. Distribution of moisture content in wet chalk cakes, 
and in cakes dried at different speeds. 


from the inner parts of the cake outwards, as 
well as a reduction of the mean water content. 
Thus, in Fig. 9, when the whizzing speed was 
increased from 980 to 1350 r.p.m., the mean water 
content of a “through”? sample decreased, but 
water content at the inner part of the cake de- 
creased more, while the outer part of the cake 
gained water. A similar effect is noticeable in 
other cases, such as at 2500r.p.m. on Figs. 10 and 13. 
The shape of drying curve is the same irres- 
pective of the depth chosen for taking the 
samples (Fig. 12) or whether it is constructed 
from the mean moisture contents of the inner, 
centre or outer layers of the cake (Fig. 13). 
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Fig. 14 is the result of tests on a chalk cake 
which was run submerged at 2570 r.p.m. before 
the dry test. The highest speed running wet is 
referred to hereafter as the formation speed. 
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Figs. 9-13. Variation of residual moisture with speed of 
basket for a chalk cake formed at 2000 r.p.m. : 





Fig. 9. Samples from the top of the basket. 


The result of using Teepol to lower the surface 
tension of the formation and wash waters is 
shown in Fig. 15, the curve’s shape remaining 
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Fig. 10. Samples from the middle of the basket. 

unchanged although the moisture contents are 
reduced. The effect of the formation speed on 
curve shape is illustrated in Fig. 16, which repre- 
sents the mean of the water content of through 





samples taken from the mid-height of cakes 
formed at different speeds. The formation speed 
abscissa coincides with a change in the effect of 
increasing the speed on the rate of cake drying 
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Samples from the bottom of the basket. 





Fig. 11. 


and hence in the curve shape. The part before 
the formation speed is steeper and concave 
downwards, whereas thereafter the gradient 
decreases, the curve being convex to the abscissa 
axis. It seems that the cake starts a consolidation 
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Fig. 12. 


at this speed, which presumably results in lower 
porosity and greater ability for the cake to retain 
its water. 

Curve B in Fig. 15 is from tests on cakes formed 
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in the small cell, spun wet at 2500 r.p.m., then 
spun dry at the prescribed speed for each sample. 
The form of curve is the same for the hydro- 

















5 


ENTAGE MOISTURE CONTENT 








“§ 


| | | 
10 





1S , 22 
SPEED, NKIO® r-p.m. 
Fig. 13. Mean of different heights at same radial positions. 


extractor cakes and small cell cakes, when spun 
dry at speeds below the formation speed. 
The inner layer of the kieselguhr cakes always 
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Fig. 14. Variation of residual moisture with speed os 
basket for a chalk cake formed at 2570 r.p.m.; samplef 
from mid-height of basket. 





rpm. 


contained less moisture than the rest of the cake. 
While the centre was of higher moisture content 
at the early stages of drying, it gave way to a 





gradual increasing moisture content distribution 
from inner to outer faces of the cake when run 
at higher speeds (Fig. 17). For comparison the 
broken curves are included, to indicate the slight 
changes with speed when running wet. Curve C 
in Fig. 17 is for a kieselguhr cake which was 
formed and washed at 1800 r.p.m., whereas the 
other curves are for a cake formed and washed 
at 2700 r.p.m. The formation speed had no 
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Fig. 15. Variation of residual moisture with speed fcr 

chalk cakes formed at 2570 r.p.m.; A: cake formed in 

hydroextractor from water. B: cakes formed in small 

cell from water. C : cake formed in hydroextractor from 
water-Teepol solution. 








effect in this case, the invariable convexity to the 
abscissa axis suggesting that the kieselguhr always 
packs in the form similar to chalk above its 
formation speed. 

When a starch cake was washed at 2500 r.p.m., 
then tested dry, the results in Fig. 18 were 
obtained. The change in the water content was 
insignificant up to a speed of about 2000 r.p.m. 
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where the cake began to dry rapidly, especially 
at the inner layer. When the surface tension of 
formation and wash water was halved by Teepol 
additions, the cake showed a similar drying 
behaviour (Fig. 19). 
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Fig. 16. Effect of formation speed on the variation of 

moisture content with speed of basket for cakes formed 
from different masses of chalk. 
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Fig. 17. Variation of moisture content with speed for 
kieselguhr cakes, submerged and dry. 


The simplest consideration of the effect of 
rotational speed on residual moisture content in 
a whizzed cake may be based on the balance 





between centrifugal force and surface tension on 
the moisture retained on the surface of the packed 
material. Consider a mass V (1 — ¢) o of material 
of specific surface S and packed volume V spun 
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Fig. 18. Variation of moisture content with speed for a 
starch cake formed at 2500 r.p.m. 
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Fig. 19. Effect of formation speed and surface tension of 
formation solution on starch cakes. 
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at a mean effective radius r and angular velocity 
w. The mean perimeter of the wetted surface in 
the bed is [S V (1 — «) o/L], if it is assumed that 
the surface can be characterised by a perimeter 
and effective length L. Then if the mass of 
liquid retained on the surface of the solid is m, 
the force balance gives 


y. SV (1 —e)e 
rw L 


The liquid to solid ratio is given as 


m _ 
Vil—ej)eo 


rw* L 


ys 


4n? Ln? r 


or R (2) 
If the solid absorbs water to give a residual 
content R., equivalent to the items (a) and (c) 
mentioned earlier, which cannot be extracted by 
centrifugal force ; then the measured ratio R is 
given by 
ys 

$n? rn? 


R= R,4 


(2a) 


This relation suggests that a plot of (R — R,) 
against 1/n* should be a straight line passing 
through the origin of coordinates, if R, is con- 
stant. This is not found to be so, the R data for 
chalk failing to conform to a simple power 
function of n, and the data for kieselguhr con- 
forming to Ra (1,/n)¥, where y = 1-0 (Fig. 20). 
In Fig. 20 data for chalk, kieselguhr and starch 
cakes formed at 2600 r.p.m. have been plotted 
as log (R — R,) against log n’, the value R, 
being varied from zero to R,, the residual 
moisture content for cakes spun dry at 2600 r.p.m. 
If R, = 0, the plot should give a straight line 
from equation (2a) of gradient — 2.0. As R, is 
increased, the gradient of the plot increases. A 
straight line of gradient — 2-0 is not obtained 
within the range 0 < R, < R,,, probably due to 
the variation of L, R, with n. Within the present 
experimental ranges of speed and cake geometry, 
it cannot be assumed that either R,, can be taken 
as R., or that any nominal value of R, is likely 
to reduce the data to the desired straight line. 
The form of the curves in Fig. 20 suggest that 
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the variations of R, and effective surface S vary 
sufficiently with compressive stress on the cake 
to render the equation (2a) valueless as a corre- 
lation function. The single effective length L is 
probably inadequate and should be replaced by a 
distribution function describing the effective mass 
and position of removable moisture. 

It is possible to devise linear plots which may 
be of value in characterising a material. When 
the experimental results were plotted as log 
(R — R,,) against log (n,,° — n’) for data at n’ 
less than the formation speed, straight lines were 
obtained in the cases of chalk and kieselguhr 
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Fig. 20. Logarithmic plot of residual moisture content 


against speed for different R, values varying from zero 
to R,,. 


(Fig. 21). The equation representing these lines 
is 


(R — R,) =C.(n,, —n)* (3) 


whence 


when 2 <1, the negative gradient dR/dn will 
increase with increasing n, as with chalk which 
has an R, n curve convex upwards (Figs. 14 and 
15) for n below the formation speed. When 
a2<1; the negative gradient dR/dn decreases 
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with increasing n, which is the case with kiesel- 
guhr, the R, n curve being concave upwards with 
respect to the abscissa axis. Thus the same type 
of equation (3) can represent the cases of both 
chalk and kieselguhr cakes by straight lines in 
Fig. 21. 
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Fig. 21. Logarithmic representation of equation (3). 


Discussion o¢ RESULTS 


The tests described above lead to the following 
conclusions : 

It is unlikely that a general equation can be 
obtained to cover the changes in porosity occurr- 
ing in a cake during its initial compression when 
submerged, as abrupt packing changes seem to 
occur. 

It is unlikely that a general function can be 
derived for materials of the type studied to 
describe the de-watering of cakes by whizzing. 
When cakes are whizzed above their formation 
speed, defined as the highest speed attained when 
submerged, the drying curve differs from that for 
lower speeds, probably due to consolidation. The 
drying curves for specific materials running at 





speeds below the formation speed may be repre- 
sented by empirical equations such as equation 
(3) which lead to linear plots of possible value. 
The speed range below the formation speed is 
the probable range of interest. 

The apparent compressive effects and packing 
distributions make it impossible to use a simple 
surface tension function for correlating de- 
watering data, even when using data from mean 
positions in hydroextractor cakes. 

The porosity distributions considered propor- 
tional to the moisture distribution in submerged 
cakes pass through maxima in a radial direction, 
suggesting that the permeability thus does not 
follow a smoothly decreasing porosity with 
increasing radius. Though the methods of cake 
formation in small cells and hydroextractor 
baskets show little influence on the mean per- 
meabilities [2], the porosity distributions may 
differ. The data in Fig. 2 show that formation 
mechanism has no effect in the small cell, differing 
from the results in the hydroextractor. The 
latter are probably dependent on whether the 
cake is formed by a multi-layer technique. 
Though the permeability differs in vertical and 
horizontal position even in cakes of fine particu- 
late material, it seems adequate to assess cakes 
in terms of the mean values in practice. 


NOMENCLATURE 

L = effective length of packed sample em 

m = mass of water in sample K 

n = speed of revolutioa ; n,, maximum test 
speed ; (rev. per sec.) ; n’ 

R = moisture content taken as ratio of water 
to solid in cake ; R,, minimum content 
attained at highest test speed ; R, resi- 
dual content. 


rev. per min. 


r = radius from axis of rotation em 
S = effective specific surface of packed 

material sq. em/g 
V = volume of packed sample ce 


y = surface tension dyne/em 
«€ = porosity of packed bed 

@ = density of solid g/ce 
w = rotational speed radians /sec. 


Il 
| 
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Le réacteur étagé—outil fondamental 
de Génie Chimique 


V. R. Lecvierc* 
Institut du Génie Chimique, Université de Toulouse 


( Recewwed 29 July 1953) 


Summary—The continuous flow stirred tank reactor (C.F.S.T.R.) is especially suited for the 
experimental study of reaction kinetics in homogeneous systems; the steady conditions 
achieved in the reactor are connected directly to the laws of chemical change and not through a 
differential equation. Experimental determination of the steady state can be performed with 
increased accuracy. 

In many cases, kinetics in a C.F.S.T.R. remain independent of the geometrical shape ; the 
results can be readily extrapolated for industrial design. A convenient graphical method is 
proposed for the accurate determination of the optimum volume ratio between the various 
compartments. While the conditions at the start or the stop of the reactor remain very difficult 
to establish by theory, with a laboratory model their determination by direct experimentation 
is very easy. 

A convenient design of laboratory set-up is proposed. 
a type of equipment will be very useful for the development of chemical engineering kinetics. 


Preliminary results confirm that such 


Résumé— Développant la théorie du réacteur étagé continu de Densicn, lauteur insiste sur 
les avantages de cet appareil pour l'étude expérimentale de la cinétique d'une réaction homogéne : 
l'état stationnaire réalisé dans chaque étage reste une fonction directe au lieu d'une fonction 
différentielle de la loi cinétique. La détermination expérimentale de cet état stationnaire peut 
étre rendue trés précise. 

Dans bien des cas, les résultats du réacteur étagé sont indépendants de la forme géométrique : 
ils sont alors directement extrapolables a l’'application industrielle. L’auteur présente une 
méthode graphique commode et précise pour le calcul du volume optimum des divers étages. 
Tandis que l'établissement théorique de la mise en régime ou de l’arrét reste toujours trés pénible, 
lessai au Laboratoire permet de préciser aisément toutes les circonstances accompagnant les 


états transitoires. 


L’auteur propose finalement un type simplifié de réacteur de laboratoire dont les résultats 
justifient lintérét de ce type d’appareil pour les études de Génie Chimique. 


Le calcul des appareils industriels peut étre con- 
sidérablement simplifié si on étend la notion 
d’opération fondamentale aux appareils eux- 
mémes, plutoét qu’A la nature chimique de la 
réaction qui s'y produit. Mais ce calcul exige 
toujours un assez grand nombre de données 
thermodynamiques et cinétiques ; la découverte 
et l'utilisation industrielle constante de nouveaux 
produits font que ces données sont trés fréquem- 
ment manquantes. Les méthodes d’approximation 
sont d’un certain secours, mais ne sont jamais 
qu'un pis aller, et sont inapplicables ou inexis- 
tantes dans bien des cas (surtout en cinétique 
chimique). 


* Tué accidentellement le 23 Septembre, le Professeur V 
que nous publions. 


Le probléme se pose donc, et de facon sérieuse, 
de savoir comment I’Ingénieur va se procurer les 
données indispensables au calcul des processus 
d’échanges intra-moléculaires. 

Nous nous proposons ici d’examiner un type 
d’outil qui nous a paru spécialement propre a 
certaines études de Génie Chimique, de montrer en 
quoi il est spécialement apte a fournir des données 
nécessaires A toute une série de calculs industriels, 
et nous proposerons un modéle type d’appareil 
peu coiteux et d’une précision suffisante. 

La théorie du réacteur étagé (en anglais : Con- 
tinuous flow stirred tank reactor) a été établie 
par DenBiGu [1]. 





. R. Lecierc n‘a pu revoir lui-méme les épreuves de l'article 


Les éditeurs déplorent cet accident qui prive l'Université Frangaise d'un chercheur passionné des 


problémes du Genie Curmigur ; ils s’excusent des erreurs que la mort de l’auteur n’a pas permis de corriger dans le 


texte ci-dessous. 
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V. R. Leciterc: Le réacteur étagé 


Rappelons le principe de l'appareil. Le réacteur 
étagé est constitué par un certain nombre d’étages 
(Fig. 1) alimentés en continu, et débitant en con- 


IF 




















Fig. 1. 


tinu, & débit constant et égal pour chaque étage. 
Chaque étage est isotherme et muni d’un puissant 
agitateur. Si nous laissons de cété le cas des 
réactions a pression constante et volume variable, 
les conditions qui suivent sont nécessaires et 
suffisantes pour que chaque étage puisse étre 
caractérisé chimiquement par un régime per- 
manent et une composition homogéne (qui est 
la composition a la sortie) : 


(a) Constance des débits d’entrée et de sortie, 

(b) Constance du volume occupé par le fluide 
dans létage (égalité des entrées et des 

sorties), 

Agitation suffisante pour que le temps de 

petit le de 


(c) 


mélange soit devant temps 
passage, 


(d) Constance et uniformité de la température. 


Si, par ailleurs, la réaction n'est pas influencée 
par les parois, il est évident que toutes les con- 
ditions ci-dessus peuvent étre réalisées independ- 
amment de la forme ou des dimensions de I '¢tage. 

Cette caractéristique du réacteur ctagé fait 
qu'il constitue un appareil de choix pour la 
détermination des constantes cinétiques et, en 
méme temps, un outil fondamental d'études de 
Génie Chimique, particuliérement apte a faciliter 


outil fondamental de Génie Chimique 


le passage direct de l'étude de laboratoire a la 


realisation industrielle. Nous allons examiner 


ces deux aspects. 


DETERMINATION DES CONSTANTES CINETIQUES 


L’étude d'une réaction en solution exige des 
déterminations : 

(a) De lordre de réaction (a l’origine et en fone- 
tion de l'état d’avancement), 

De la constante de vitesse, 

De linfluence de la constante diélectrique du 


solvant et de la force ionique. 


(b) 
(¢) 


Ces déterminations se raménent finalement a 
des mesures de composition chimique et de 
temps. 

La plupart des données de la litterature pro- 
viennent de déterminations faites avec un réacteur 
simple (batch process). Lorsqu’il s‘agit dune 
réaction en solution, on opére généralement sur 
un volume de liquide assez grand, obtenu par 
mélange rapide de deux réactifs. 

Pour suivre la composition en fonction du 
temps, on effectue, & des intervalles de temps 
choisis, des prélévements qu’on analyse (sauf le 
cas relativement rare olf on peut suivre en con- 
tinu la composition par une méthode physique). 
Ces prélévements sont fréquemment effectués a 
l'aide d'une pipette. II résulte de ce mode 
opératoire une incertitude assez grande tant sur 
lorigine des temps que sur la “date” des 
prélévements. De plus, linterprétation des 
résultats sera assez difficile la réaction 
d’ordre fractionnaire et si la force ionique varie 


sl est 
au cours de la réaction (5-16). 

Nous voulons montrer que le réacteur étagé offre 
sur le réacteur simple des avantages notables quant 
d la précision et d la simplicité des déterminations. 

Considérons une réaction (irréversible pour 
simplifier) 
>~R 


nA+mRPR 


d’ordre « par rapport a A et d’ordre 8 par rapport 
a B c’est-a-dire dont la vitésse obéit a la loi : 
da 
= — k a® LP 
dé 


ou a et b représentent les concentrations instan- 





tanées en A et B. 
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I] s’établira dans chaque étage un état station- 
naire correspondant & des concentrations a, a, . . . 
a, et b, b,...b,. L’alimentation du réacteur est 
caractérisée par des concentrations a, et by. Soit 
F le débit de cette alimentation, V,, V,...V,, 
les volumes des étages. 

Il entre dans le ler étage, par unité de temps, 
a,F moles de A et il en sort, dans le méme temps, 
a,F. La différence est F(a, —a,) pour le 
volume V’,, et par unité de volume F'/V, (a,—a,) 
pour le ler étage. Pour le néme étage, cette 
différence devient F/V, (a,_, — 4,). 

La différence entrée — sortie est égale au nombre 
de moles de A disparues par unité de temps et 
unité de volume, du fait de la réaction. Pour 
l'étage 1 de concentrations stationnaires a, et },, 
elle est par définition : 

da 


Sa k af bf 


D’ou les relations, pour l’étage 1 : 
— kg? BP 

pr (dg —a,) = kay by 
1 


pour l’étage n (si l’ordre ne varie pas) : 


n-1 a,,) =k a, by 

Ce type d’équations contient des grandeurs a 
déterminer expérimentalement et les inconnues 
k, a, B. Ceci permet, au prix d’un nombre con- 
venable d’expériences de déterminer l'ordre 
(méme fractionnaire), et la constante de vitesse. 
DenbiGu [1] a fait remarquer l’intérét de ce type 
d’équation : il n’y a pas a effectuer d’intégration, 
avantage considérable pour les réactions com- 
plexes. 

Par ailleurs, la précision du maintien de l'état 
stationnaire peut étre trés élevée (meilleure que 
0-1%) car elle ne dépend que du réglage des 
débits et d'une agitation suffisante; tout se 
passe comme si la réaction était “ gelée”’ a un 
état donné. La précision des mesures de con- 
centration peut également étre trés élevée 
puisqu’il est possible d’effectuer autant de préle- 
vements que l’on veut pour définir l'état per- 
manent—on peut done faire la moyenne de 
nombreuses analyses chimiques. D’autre part, le 
temps correspondant au prélévement est ici le 





outil fondamental de Génie Chimique 


rapport V'/F du volume occupé au débit de 
lalimentation. Ce temps peut étre déterminé 
statistiquement et non par une seule mesure comme 
dans un réacteur simple. 

Enfin, en recevant le liquide sortant de l’appa- 
reil dans un récipient contenant un réactif con- 
venable, il devient possible d’arréter instantané- 
ment la réaction. D’autre part, le réacteur étagé 
permet de choisir, et de maintenir constantes, 
diverses conditions opératoires qui influent sur la 
vitesse de réaction. DenBicn [1] écrit a ce sujet : 


“La réaction peut étre ¢tudiée dans des 
conditions rigoureusement constantes de con- 
centration en catalyseur (ou en radicaux 
libres) ou de force ionique, et ceci & n’importe 
quel stade de la réaction.” 


Précisons que la force ionique peut étre main- 
tenue constante non seulement dans chaque 
étage considéré individuellement, mais dans 
l'ensemble des étages. II suffit pour cela d’effectuer 
des additions convenables d’ions & chaque étage. 
On peut ainsi, du début a la fin, étudier une réaction 
a force ionique constante, et choisie, ce qui est 
pratiquement impossible avec un réacteur simple. 

Il semble, & premiére vue, que nous venons 
d’énumérér des avantages qui appartiennent non 
au réacteur étagé mais aux réacteurs continus en 
général; c’est-d-dire qu’on pourrait aussi bien 
utiliser un réacteur tubulaire pour cette étude 
cinétique. Mais le réacteur tubulaire présente des 
inconvénients trés graves. En effet, il est facile 
de montrer que pour un tel type de réacteur, 
fonctionnant en régime laminaire, chaque tranche 
mince est caractérisée par un état stationnaire 
mais cet état stationnaire n'est pas homogéne. (Le 
réacteur tubulaire est par contre un instrument 
intéressant pour déterminer commodément le 
degré d’influence des parois.) 


Proget D’APPARILLAGE INDUSTRIEL 


Dans ce qui précéde, nous avons montré que le 
réacteur étagé, & un ou plusieurs étages, est un 
outil de choix pour la détermination des con- 


stantes cinétiques. II reste encore l’outil de 
choix pour les études de Génie Chimique destinées 
a permettre par extrapolation la réalisation indus- 
trielle d’une réaction étudiée au laboratoire. 
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Dans le cas d’une réaction effectuée dans un 
réacteur étagé, le probléme immédiat est la 
détermination du volume des étages et de leur 
nombre. De plus, il faut prévoir les conditions 
de mise en régime et d’arrét de l’appareillage. Nous 
avons déja insisté sur le fait que toutes les condi- 
tions caractérisant le fonctionnement du réacteur a 
Gtages sont indépendantes de la forme ou des dimen- 
sions des étages. L’appareil de laboratoire est 
directement transposable aux dimensions exigées 
par le tonnage de fabrication envisagé. 


Calcul des dimensions du réacteur 


E_prinGE et Piret [4] ont donné une méthode 
graphique qui permet de calculer le nombre des 
étages nécessaires & réaliser une réaction donnée 
(V /F étant fixé) ou de calculer le rapport V/F 
si on se fixé le nombre des étages. Mais pour des 
raisons de commodité, il se sont limités au cas du 
réacteur A n étages identiques. 

L’intérét qu'il y a A faire varier le volume 
relatif des étages d’un réacteur industriel n’ayant 
pas, & notre avis, été suffisamment mis en lumiére, 
nous commencerons par montrer cet intérét, 
avant d indiquer les ressources du réacteur étagé de 
laboratoire pour le calcul de l'appareil industriel. 

Il est possible d’obtenir un gain de rendement 
si on accroit le volume des étages & mesure que la 
réaction progresse (ce qui revient A accorder 
davantage de temps A une réaction qui se 
ralentit). 

Soit une réaction irréversible A+ B-+>R 
d’ordre 2 avec concentrations initiales égales 
a, = b,. La méthode du bilan matériel indiquée 
précédemment permet d’établir les équations des 
deux premiers étages : 


kaha, +a, — ay =0 (1) 


V. 
ka,’ +a, — a, =0 


(2) 
les étages suivants étant également caractérisés 
par d’autres équations du 2éme degré, la con- 
stante de chacune d’elles étant la racine de la 
précédente. 

Denpicnu [la] a fait remarquer que si on 
néglige le terme a, de l’équation (1) et le terme 





a, de l’équation (2) il est facile de tirer une 
valeur de a, qui est : 


= F, a - 
le (p VV 


expression qui passe par un minimum (corres- 
pondant a un maximum de rendement) si V7. =2V, 
~ Pour une réaction d’ordre 3, on aurait V,=3 V, 
~— En premiére approximation le rapport optimum 
des volumes serait celui de l'ordre de la réaction. 
Nous avons cherché a établir une méthode 
commode et précise de détermination de cet 
optimum et nous nous sommes arrétés a une 
solution graphique. L’équation 2 peut s’écrire : 
‘4 


a, =k 2a, +a, (3) 


et l’équation (n) serait : 
(4) 


Ces équations font apparaitre une relation 
linéaire entre le volume V d’un étage considéré 
et la concentration a laquelle on l’alimente (F 
étant fixe) si nous nous imposons la concentration 
de sortie. 

Il est alors possible de tracer une famille de 
droites divergentes qui exprimeront des relations 
linéaires entre V'/F et la concentration d’alimen- 
tation ceci pour diverses concentrations fixes de 
sortie. 

Supposons que nous voulions pour la réaction 
précédente réaliser la transformation de A de la 
concentration 0-7 & la concentration 0-5 en 2 
étages. Considérons la Fig. 2. La droite D, 
représente |’équation : 


a, =k "20.5 + 0-5 


et le point I correspond a une valeur de V/F qui 
impose une concentration d’alimentation 0-52. 
Si nous menons par ce point une paralléle & 
la droite D, (a, = k V,/F.0.52* + 0-52) nous 
atteignons l’ordonnée 0-7 au point I’ dont l’ab- 
scisse donne pour V’/F la valeur correspondant 
au temps total de séjour pour réaliser la conversion 
de 0-7 8 0.5, Si nous recommencons |’opération en 
nous arrétant au point Ii d’abscisse 0.54 nous 
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tracerons une paralléle & la droite D, (a, = K 


nous atteindrons l’ordon- 


J 0.542 + 0-54) et 


née 0-7 en un point II’ dont l’abscisse est plus 
faible que celle du point I’. Si nous continuons 
a faire croitre les temps de passage alloués au 
dernier étage, nous verrons que le point final 
passe par une abscisse minimum, ce qui fixe le 
rapport optimum des 2 étages. 

Si, par convention, nous prenons pour unité 
les concentrations initiales a, et b,, nous obtien- 
drons, en cing essais, les valeurs suivantes pour 
les temps de passage relatifs & travers deux 
étages successifs : 





Temps de passage 
1* Aage 
(V/F), 

447 
358 
315 
245 
243 


total 
(V/P), 
687 
678 
675 
675 
683 


2° dlage 

(V/F)y 
240 
820 





Le tracé de la courbe (V/F), en fonction de 
(V /F), donne pour valeur minimum (V /F),— 674 
avec (V/F), = 380 et (V/F), = 294 ce qui 
correspond a un rapport optimum d’environ 1-3. 

La méme opération effectuée avec une réaction 
d’ordre 3 donne 1-7 environ pour le rapport 
optimum des valeurs de sorte que le rapport 
optimum du volume des étages varie A peu prés 
comme la racine carrée de l’ordre de la réaction. 

Il est facile de vérifier la position du minimum 
en comparant un réacteur A étages égaux, un 
réacteur A étages inégaux et un réacteur tubu- 
laire.I] n’y a qu’a appliquer les formules classiques 
pour le réacteur tubulaire, et la série d’équations 
du 2éme degré pour n étages égaux ou inégaux. 

La Figure 3 donne un exemple de cette com- 
paraison. 

Cette méthode est évidemment générale,pourvu 
que l’on connaisse l’ordre et la constante de 
vitesse de la réaction ; il reste encore possible de 
l’'appliquer si l’ordre et la constante de vitesse 
varient au cours de la réaction, soit de par la 
nature de la réaction, soit & cause d’additions 
variées. 





On peut choisir des conditions opératoires aussi 
complexes que l’on veut (choix dicté par des 
études cinétiques préalables). L’usage d’un 
réacteur étage de laboratoire permettra toujours 
de déterminer, pour une zone peu étendué de 
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Fig. 2. Schéma de la méthode graphique pour calculer le 
rapport optimum des volumes des étages. 


conditions opératoires, la loi de variation des 
coefficients cinétiques; il permet évidemment 
de déterminer graphiquement le rapport de 
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Fig. 3. Comparison du reacteur tubulaire ————. 
Comparaison du reacteur étagé (étages egaux) — — —. 
Comparaison du reacteur etagé (étages inegaux 
Reaction A +B>R. ay =by. k = 107%. 


volumes convenant a ces conditions complexes 
et de vérifier, au laboratoire, ce que donneront 
les appareils industriels. 
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Une fois établi & l'aide des données de labora- 
toire, le graphique permet de choisir le nombre 
d’étages nécessaires; ce choix dépendra des 
contraintes que l’on s’imposera par des considéra- 
tions d’ordre économique. Dans beaucoup de 
cas, partant d'une concentration initiale donnée, 
il suffit d’aboutir a une concenfration finale, plus 
ou moins rigoureusement imposée : le graphique 
permet la détermination du nombre d’étages a 
choisir; il montre immédiatement les avantages 
que présentent des étages inégaux pour contre- 


balancer les complications dans la construction. 


Mise en marche et arrét 


Pour mettre un réacteur étagé en marche, on peut 
utiliser procédés divers : 


(a) Tous les étages sont remplis par un fluide inerte 
ou par le solvant pur, 

(b) On remplit chaque étage avec l’alimentation 
initiale, 

(c) On remplit chaque étage en y mélangeant des 
solutions titrées de fagon & produire la com- 
position prévue pour l’état stationnaire. 


De toutes fagons, la prévision de l’atteinte de 
l'état stationnaire (nécessaire pour savoir ce qui 
sortira de l'appreil pendant ce temps) présente 
trés vite de grosses difficultés mathématiques 
(voir & ce sujet Prret et Mason [2)). 

Les qualités d’extrapolation du réacteur de 
laboratoire permettent une prévision compléte 
des compositions chimiques pendant la mise en 
marche et l’arrét, dans des conditions en principe 
aussi compliquées que l'on veut (toujours bien 
entendu si la paroi de l'appareil est inactive 
chimiquement). 


Reactreur EraGe pe LABORATOIRE 


Les considérations précédentes nous ont amené 
a considérer comme un probléme trés intéressant 
celui de la réalisation d'un bon réacteur de 
laboratoire. Divers modéles ont déja été proposés 
en particulier celui de Dennicu [3] et celui de 
Piret et Exvprince [4]. Ces modéles ont pour 
caractéristique d’étre entiérement clos pour per- 
mettre une connaissance exacte du volume de 
étage. 

Nous avons cherché a établir un modéle ouvert, 





peu cofiteux a réaliser, facile & manipuler et de 
precision suffisante. Mais pour justifier les 
considérations énoncées plus haut nous avons 
choisi parmi plusieurs modéles possibles celui qui 
avait la forme la plus différente des modeéles déja 
décrits : le modéle retenu est en fait un tube 
de Thiele agrandi. 

La Figure 4 donne les caractéristiques essen- 
tielles de sa construction. Si la tige de lagitateur 








Fig. 4. Corps en pyrex 30/34 mm. 
Alimentation par vases de Mariotte. 
Chauffage par rayonnement. 

Thermo régulateur & mercure-Soudeé. 

Heélice 4 pales, arbre ¢ = 4mm. 

Repére gradué ()" = 340 cc). 

Robinet de deversement dans le 2éme étage. 
Pouite capillaire pour prélévements en cas de 
fonctionnement a plusieurs étages. 


est assez fine, il n'y a pas formation de vortex 
appréciable et on peut s’assurer de la constance 
du volume comme pour une fiole jaugée. Le 
fonctionnement de l'appareil a été soumis aux 
vérifications expérimentales suivantes : 


(a) Agitation—La rapidité d’homogénéisation 
a été vérifiée par remplissage de l’étage avec une 
solution de soude caustique N/5 et addition dune 
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goutte de phénol phtaléine, (on peut utiliser toute 
réaction colorée instantanée) et par addition a 
de l'eau distillée de bleu de méthyléne concentré. 
Le temps de mélange est de 3 secondes + } 
seconde, ce qui est suffisamment petit pour 
utiliser des temps de passage de l’ordre de 300 
secondes et plus. 
(b) Débits—La constance des débits des vases de 
Mariotte et du réacteur a été vérifiée par remplis- 
sage chronométré de fioles jaugées. 
Le débit en cc./minute est d’une bonne repro- 
ductibilité (écart type de 0.2%). 





| 
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0-002 








0-00! : ? ; 
TEMPS EN HEURES 














Fig. 5. Atteinte de l'état stationnaire. 
[NaOH], = 0-00498 moles/litre. (CH,CO,C,H,)] = 0-225 
Temps de passage = 22 min. 


Pour les premiers essais effectués & 25°C des 
moyens trés simples suffisant a régler la tem- 
pérature A + 0-1°C, Ces premiers essais regardés 
comme satisfaisants, nous avons essayé l'appareil 
en y réalisant une réaction connue. 





Nous avons choisi la saponification de l’acétate 
d’éthyle en solution diluée, réaction dont la 
constante de vitesse a été établie par des travaux 
classiques, et chose intéressante, vérifiée par 
DENBIGH [3] avec un modéle fermé de réacteur 
étage. 

Les analyses ont été faites en pesant des 
quantités donées du produit sortant, recu dans 
de l’acide chlorhydrique titré N/100 et titrant 
l'excés d’acide en présence d'indicateur de 
neutralité. 

La Figure 5 donne la courbe d’atteinte de 
état stationnaire. Ceci permet de calculer la 
valeur de la constante de vitesse k. On trouve 
10-8 10? + 0-2: Denbicnu a précisément trouvé 
cette valeur, qui correspond également aux 
travaux classiques de WALKER et REICHERT. 

Le temps d’atteinte de l'état stationnaire est 
aussi pratiquement celui de DENBIGH, environ 8 
fois le temps de passage (le degré d’approximation 
est fonction du temps). 

I] est remarquable que l'appareil de DENBIGH 
et le nétre, de forme complétement différente, 
donnent le méme résultat pour la méme réaction. 


L’appreil de Piret qui est aussi différent de l'un 
que de l'autre donne les résultats prévus d’aprés 
les mesures statiques pour une autre réaction. Ceci 
montre que le réacteur étagé est bien un appareil 
extrapolable pourvu que les conditions que nous 
avons énoncées soient respectées. 


Les opérations expérimentales citées ici ont 
été en grande partie effectuées par Mr. Denys 
Vesco A qui nous exprimons ici nos sincéres 
remerciements. 
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Etudes d'equilibres liquide-vapour 2 


(a) Systeme Benzéne-Diméthylformamide sous un atmosphére absolu 
(6) Systéme Méthanol-Diméthylformamide sous un atmosphére absolu 


A. DELZENNE 


Laboratoire de Recherches des Etablissements Kuhlmann, La Madeleine, Paris 


(Received 1 Sept. 1953) 


Summary— Experimental study of the liquid-vapour relationships in two binary systems, with 
quite different boiling points. Comparison with Gibbs-Duhem theoretical equation. 


Two different experimenjal methods have been necessary to get reliable results whose validity 


is checked by various criteria. 


No previous figures have been found in the literature for these binaries. 


Résumé—Détermination expérimentale du diagramme d'équilibre liquide-vapeur des deux 


binaires : 


Benzéne-Diméthylformamide 
~ Méthanol-Diméthylformamide 


Les résultats sont comparés avee ceux déduits de l'équation théorique de Gibbs-Duhem. 


L‘intérét de ce travail réside dans l'étude de deux systémes & points d’ébullition éloignés par 
deux méthodes différentes et dans la comparaison entre différents critéres de verification. 


Les deux binaires étudiés ici, n'ont fait lobjet, 4 notre connaissance, d’aucune publication. 


1. Purete pes ConstIruANTs 
(a) Benzéne 
Le produit de départ est un benzéne technique 
purifié par des procédés chimiques puis par 
rectification soignée. 

Le processus de purification a été le suivant : 

(a) lavage du produit technique par l’acide 
sulfurique concentré jusqu’A ce que l'acide ne 
soit plus coloré. 

(b) lavage par une solution de carbonate de 
soude puis & l'eau pure. 

(ec) séchage par SO*Na? anhydre. 

(d) refroidissement de facon A ce que les } 
environ de la quantité traitée soient solidifiés, le 
dernier quart étant rejeté. Cette opération est 
répétée trois fois; on élimine ainsi les traces de 
toluéne. 

(e) nouveau séchage, puis rectification sur 
sodium dans une colonne A haut pouvoir sépara- 
teur. 

On arrive de cette facon A un produit ayant un 
degré de pureté trés élevé comme I'indique le 
tableau 1. 





(b) Méthanol 

Le méthanol est obtenu pur par rectification dans 
une colonne A haut pouvoir séparateur. Dans le 
Tableau 1, on compare les caractéristiques du 
produit utilisé avec celles données par littérature. 


(c) Diméthylformamide 

La diméthylformamide est obtenue a partir du 
formiate de méthyle et de la diméthylamine 
anhydre suivant la réaction : 


{| 
CH-O-CH, + (CH,)*NH +CH,OH +CH-N(CH,)? 


Les deux constituants du mélange réactionnel 
(méthanol et diméthylformamide) sont séparés par 
distillation. La diméthylformamide est ensuite 
épurée par rectification sous un vide de 15 mm 
Hg. ce qui correspond A une température de 
passage de 48°. 

Les constantes physiques des produits obtenus 
aprés les traitements d’épuration figurent dans 
le tableau 1. 

Les indices de réfraction ont été déterminés au 
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Tableau 1 





Benzéne 


Experi mental 





Indice de réfraction np*° 


(1) 


Littérature 


1:50142 


Méthanol Diméthylformamide 





Littérature Littérature 


Expérimental 


Expérimental 





1-:33118 
[1] a 14°S 


1-4269 
[2] a 25° 





Densité d,?° 0-879% 


(1) 














0-7928 0-9445 


[2] & 25° 














moyen du réfractométre Fery-VAURABOURG don- 
nant la 4éme décimale exacte. 

Les densités ont été déterminées au moyen d'un 
picnomeétre. 


2. Tensions DE VAPEUR 
Le Tableau 2 reproduit les valeurs des tensions 
de vapeur des constituants purs trouvés dans la 
littérature, tout au moins en ce qui concerne le 
benzéne et le méthanol. 

Par contre, nous n’avons trouvé dans la littéra- 
ture que bien peu de données relatives A la 
diméthylformamide. Seuls sont mentionnées les 
points d’ébullition sous 758 mm Hg et 39 mm Hg 
[3] et il a été jugé préférable de déterminer 





expérimentalement la courbe de tension de 
vapeur de la diméthylformamide. 

On a utilisé pour cette détermination un 
ébulliométre de Swietolawski, la partie supérieure 
du condenseur étant reliée d’une part a un mano- 
métre A mercure précis, et d’autre part a une 
pompe & vide. Entre la pompe et |’ébulliométre 
on a placé un réservoir de 10 litres servant de 
régulateur de pression. On introduit dans 
l’ébulliométre, un mélange de composition déter- 
minée et quand la pression désirée est atteinte, on 
note la température indiquée par un thermométre 
de précision. 

Les résultats sont reassemblés dans le Tableau 2 
et représentés sur la Figure 1. 


Tableau 2 





Pression en mm Hg | 


aie | 10 


| 
Température en °C | — 31-7 | —196 | —115 | - 2-6 
| | 





oo | "200 


| +76 


| 
+ 15-4 <a +422) + 60-6 





Pression en mm Hg 


Température en °C — 16-2 











Méthanol (4) 


——— |} —— ——_ , —___—_. — 


40 «| 660 | «(100 


fe > 
+121 | +212 


| 200 


+ 5 


+ 34-8 | + 49-9 











— 








Diméthylformamide 





Pression en mm Hg 








- 
200 


400 760 





Température en °C + 63-5 | +78 








+105 | + 130 





+ 153-4 
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Courbe de tension de vapeur de la dimethylforma- 
mide. 


Fig. 1. 


Pour I’interpolation nécessaire a tous les calculs, 
on a utilisé la formule connue : 
B 
T 
qui se traduit par une droite sur un diagramme 
a ordonnées logarithmiques. 


log P = A — 


3. ANALYSE DES ECHANTILLONS 
L’analyse des échantillons a été faite par réfracto- 
métrie, les courbes n,*° = f(z) ayant été déter- 
minées au préalable pour les deux systémes 
étudiés. Iles valeurs correspondantes sont 
reportées sur le Tableau 3 et sur les Figures 2 et 3. 


Tableau 3 
| 
Systeme 
Méthanol- Diméthylforma- 
mide 





Systeme 
Benzéne- Diméthylformamide 


Indice de 
refraction 


"Dp 


Indice de Titre 
réfraction moléculaire 
np*” en méthanol 


Titre 
moléculaire 
en benzéne 





1-4710 
1-4835 
0-4922 
15012 
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oe 
Fraction molaire en benzene 
Systéme benzene-dimethylformamide. Courbe 


Fig. 2. 
np?” = f (x). 


4. Description Des Metuopes EXPERIMENTALES 


(a) Méthode par condensation et recirculation du 
condensat 

La détermination des équilibres liquide-vapeur 

par cette méthode a été faite en utilisant l'appareil 
d’Orumer [5]. 

Une charge totale de 400 cc. est utilisée pour 

chaque essai, le régime est établi au bout de deux 
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Fig. 3. Systéme Methanol-Dimethylformamide. Courbe 


np?” = f (2). 
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A, 8. Saturateurs F Gaine thermométrique 

C Chambre d'equilibre Ge. Bain thermostatique 

R Vaporiseur. 4. Enroulement chauf fant 

D Prise d'échantillon 
oe /a phase liquide 

&. Prise d'échantillon 
de /a phase vapeur 















































| | 
a 


Appareil pour la détermination des équilibres 
Liquide. Vapeur par la méthode dynamique. 








Fig. 4. 


heures d’ébullition environ. On préléve alors un 
échantillon de la phase vapeur condensée et un 
échantillon de la phase liquide que l'on soumet 
immédiatement a l’analyse. 

La température est enregistrée au moyen d’un 
thermométre de précision A 4- 0-1°C, 


(b) Méthode dynamique 


Comme on le verra plus loin, l'emploi de l'appareil 
d’OruMeR ayant donné lieu & des mécomptes pour 
l'étude de l'un des systémes, nous avons été 
amenés a utiliser dans ce cas la méthode dynamique 
dont nous rappelons briévement le principe. 

On introduit dans une chambre d’équilibre 
portée a une certaine température et contenant 
le mélange & étudier, le constituant le plus volatil 
sous forme gazeuse. 

L’appareillage est représenté sur la Fig. 4. Il 
comprend deux saturateurs A et B et une chambre 
d’équilibre C ; ces trois récipients, reliés entre eux, 
sont plongés complétement dans un bain d’eau ou 
@huile suivant la température que l'on veut 
atteindre. La chambre d’équilibre est munie 
d'une tubulure plongeante permettant de prélever 
un échantillon de la phase liquide, d'une deuxiéme 
tubulure permettant l’évacuation des vapeurs et 





d’une gaine thermométrique pour |’enregistrement 
de la température. 

Le constituant volatil vaporisé dans le récipient 
R est envoyé dans les saturateurs A et B et dans 
la chambre d’équilibre C, ot: les vapeurs barbotent 
dans le binaire A étudier. La phase vapeur est 
condensée & la sortie de la chambre d’équilibre 
et recueillie dans un récepteur. On maintient la 
température constante dans le bain et aprés deux 
heures de barbotage, on préléve un échantillon 
de la phase liquide et un échantillon de la phase 
vapeur condensée. 

La température est lue au moyen d'un thermo- 
métre de précision (+ 0-1°C). 


CaLcuL pDES COEFFICIENTS D’ACTIVITE — 
VERIFICATION DES RESULTATS 
EXPERIMENTAUX 

(a) Calcul des coefficients d activité 

Il est connu que la valeur des données expéri- 
mentales est vérifiée dans les meilleures conditions 
par le calcul des coefficients d’activité. Ces 
coefficients s’obtiennent au moyen des relations 
suivantes : 


- Py, 
sii P, 2, 
_ Pye 
Ys — Py ay 


(1a) 


(1b) 


oll y, et y, = coefficients d’activité 


P = pression d’expérience 
x = titre moléculaire dans le liquide 
y = titre moléculaire dans la vapeur 


et ces calculs permettent de tracer les courbes 
Y = f(2,) et Ye = f (2,) 


(b) Verification 

(a) L’exactitude de ces résultats est vérifiée en 
premier lieu au moyen des équations de Van 
Laark déduites de la relation générale de Gibbs- 
Duhem. 

Les courbes résultant de l’application de ces 
équations sont ensuite tracées en pointillés pour 
comparaison avec les courbes déduites directe- 
ment des données expérimentales (en trait plein). 
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Fraction molaire en benzene 
Fig. 5. Systéme benzene-dimethylformamide. Courbes 


d'ébullition et de rosee. 


(6) On vérifie ensuite si les pentes des tangentes 
des courbes log y = f (z,) répondent a |’équation. 


dlog y, _ (=) (“BE % (4) 
dz, zr dz, 


On choisit pour cela un point sur la courbe log 
y, =f (z,) ayant pour abcisse une valeur 
déterminée z, et on mesure la pente de la tangente 
en ce point. On mesure d’autre part, la pente de 
la tangente & la courbe log y, = f(z,) en un 
point situé sur la courbe log y, = f(z,) et ayant 
pour abcisse la méme valeur z,. Le rapport de 





ces pentes doit étre égal a — ( =] 


(c) On vérifié enfin que les courbes log Ya f (2) 
Ye 


satisfont a la relation : 
1 


| toe 2 =0 (5) 
Y2 


Cette relation s’établit & partir de l’excés d’énergie 
libre par molécule de mélange, grandeur qui est 


égale a: 
AF=RT. Z[2,\n y,) (6) 


En d’autres termes, on vérifie donc |'égalité des 
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Fraction molaire en benzene dans le lguide 


Fig. 6. Systéme benzene-dimethylformamide., Courbe 


y =f (2). 
aires comprises entre la courbe log 3 =f (z,), la 
droite horizontale d’ordonnée log “1 — © et les 
Ye 


droites verticales d’abcisses z, = 0 et 2, = 1. Ces 
aires doivent étre égales, & quelques unités pour 
cent prés, de leur valeur. 


Tableau 4 


Données expérimentales de Uéquilibre liquide-vapeur du 
systéme Benzéne- Diméthylformamide 





Temp. log log 
en °C a, Y) “1 Ye ‘1 Ye 





140 0-055 | 0-355 | 1:380 | 0-935 | 0-140 |— 0-029 
1315 | 0-100 | 0-490 | 1-260 | 1-005 | 0-110 | 0-002 
119 0-200 | 0-690 | 1/160 | 1:030 | 0-075 | 0-013 
106 0-350 | 0-835 | 1:110 | 1-060 | 0-050 | 0-025 
95-2 | 0-560 | 0-925 | 1-062 | 1/120 | 0-026 | 0-052 
91 0-655 | 0-945 | 1045 | 1/185 | 0-019 | 0-074 
85 0-800 | 0-975 | 1-029 | 1-230 | 0-012 | 0-090 























6. Resuttats ExPeERIMENTAUX 
(a) Systéme Benzéne- Diméthylformamide 
La détermination de |’équilibre liquide-vapeur de 
ce systéme a été faite au moyen de l'appareil 
d’OTHMER. 
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Fig. 7. Systéme benzene-dimethylformamide. Courbes 


log y f (x). 


Les résultats expérimentaux sont rassemblés 
dans le tableau 4 et représentés sur les Figures 
5 et 6. 

La Figure 5, concerne le faisceau d’ébullition 
et de rosée; la Figure 6 traduit la relation 
y = f(x) indépendamment de la température. 

Les coefficients d’activité calculés au moyen 
des équations (la) et (1b) sont reportés sur la 


Figure 7. 
Controle et verification 


Les valeurs de A et B calculées a partir de z, = 
0-560 sont les suivantes : 


A = 0-160 
B = 0-135 


Les courbes résultant de l’application des équa- 
tions de Van Laar, A partir de ces valeurs A et 
B, ont été tracées en pointillés sur la Figure 7. 
L’accord entre ces dernieres et les courbes 
expérimentales (en trait plein) est fort satisfaisant 
si l'on tient compte, d’une part que les équations 
de Van Laar sont établies pour une série de 
mesures isothermes et, d’autre part, que |’écart 
entre les points d’ébullition est assez considér- 
able. Une telle concordance a d’ailleurs déja été 
observée par d’autres auteurs [6]. 

La relation 4 est, elle aussi, sensiblement 
vérifiée comme l’indique le tableau 5, suivant, 
tableau dans lequel M, représente le premier 
membre de la relation, M, le deuxiéme membre. 
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Fig. 8. Systéme benzene-dimethylformamide. Courbe 


log ”" = (z,). 
Y2 


Tableau 5 





ry M, M, 


0-500 
0-267 
0-165 
0-125 


0-050 


0-491 
0-289 
0-150 
O-115 
0-052 


0-100 
0-200 
0-350 
0-560 
0-800 





Enfin on a tracé sur la Figure 8, la courbe 


v1 


log “! = f (a,) & partir des données expérimentales 
Y 


2 
L’égalité des aires hachurées, réalisées & 4.3% 
prés, montre que la courbe satisfait a la relation 9. 


(b) Systéme Méthanol- Diméthylformamide 

Les premiéres séries d’expériences ont été effectuées 
avec un appareil d’'Orumer semblable au précé- 
dent; mais on a obtenu des résultats trés dis- 
persés — (Tableau 6, Figure 9). 

Les irrégularités constatées sont dues au fait 
que l’ébullition, au lieu d’étre réguliére, se faisait 
par A coups avec formation de trés grosses bulles 
de vapeur qui perturbaient |’établissement de 
léquilibre. 

Différents remédes, tels que l’adjonction de 
corps poreux ou de tubes de verre, le chauffage 
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interne par résistance électrique, ont été essayés 
sans obtenir l’ébullition réguliére désirée. 


Tableau 6 





Température | 


"1 


1-220 
1-085 
0-781 
0-710 
0-765 
0-895 
0-935 
0-952 


0-720 
0-810 
0-840 
0-932 
0-980 
0-990 
0-996 


0-250 
0-405 
0-595 
0-700 
0-820 





—~——— 





























Fraction molare en methanol 
Fig. 9. Systéme methanol-dimethylformamide. Courbe 
d’ébullition et de rosée (Methode par recirculation-Appareil 
d’Orumer). 























Temperature 
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a % 
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0 ; O4 06 
Fraction mola:re en methanol 
Fig. 10. Systéme methanol dimethylformamide. Courbes 
d’ebullition et de rosée. (Methode dynamique). 


D'aprés ce tableau il est visible que les valeurs 
de y ne sont pas acceptables. Si abstraction faite 
de la dispersion des points expérimentaux, on 
trace une courbe moyenne pour y, et pour y,, ces 
courbes passeront dans l'un et l'autre cas par un 
minimum. Ce qui n’est pas justifié au point de 
vue thermodynamique. 

Devant cet insuccés, on a donc décidé de faire 
l'étude de l'équilibre liquide-vapeur de ce 
systéme par la méthode dynamique qui a été 
décrite plus haut. 

Les résultats expérimentaux ont été rassemblés 
dans le tableau 7, et les courbes qui en sont 


Tableau 7 
Données expérimentales de Uéquilibre Liquide-Vapeur du systéme Méthanol- Diméthylformamide 





Température en °C | 


0-300 
0-525 
0-675 
0-790 
0-875 
0-930 
0-968 
0-9924 











log %1 log Ye 





— 0-208 
— 0-172 
— 0-121 
— 0-100 
— 0-072 
— 0-041 
— 0-019 
— 0-018 
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F Tableau 8 














S 
2 





0-145 
0-210 
0-300 
0-415 
0-575 





) 
yn 











Ss 
= 





faisant entre les courbes expérimentales (en trait 
plein) et les courbes déduites des équations de 
Van Laar (en pointillés, apartir de ces valeurs de 
A et B. 

a3 34 ry, D’autre part, le relation qui existe entre les 
Fraction molaire en methanol dans le /iquide pentes des tangentes aux courbes log y, = f (z,) 


Fig. 11. Systéme methanol-dimethylformamide. Courbe | &t log vy, = f (z,) est approximativement vérifiée 
d’équilibre liquide-vapeur. comme I’indique le tableau 8. 
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Fraction molare en methanol dans /e liquide 
Fig. 12. Systéme methanol-dimethylformamide. Courbes 
log y = f (2). 


déduites, tracées sur les Figures 10 et 11. Les En dernier lieu, on a tracé sur la Figure 13, la 
coefficients d’activité, calculés comme précédem- 
ment, sont reportés sur la Figure 12. 


courbe log ”! = f(a,). L’égalité des aires hachu- 
Y 


rées est vérifiée & 3% prés. Les données expéri- 
mentales satisfont donc, de facon acceptable a 
la relation : 


Vérification et controle 
Les valeurs de A et B, calculées & partir de 7, = 
0-415, sont les suivantes: 
A= — 0-196 
B= — 0-125 
La Figure 12 montre qu’il y a un accord satis- 
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pas des systémes idéaux. 


Si l'étude du premier systéme est réalisable par 
la méthode de recirculation au moyen d’un 
appareil d’'OruMmer courant, celle du second s'est 
révélée irréalisable dans ces conditions. La 
méthode dynamique a permis de surmonter cette 
difficulté. 

Les résultats expérimentaux sur |'équilibre 
liquide-vapeur de ces deux systémes ont été 
soumis a différents critéres dérivés de la thermo- 
dynamique des solutions. Les courbes expéri- 
mentales obtenues montrent vis-A-vis de celles 
cee a ia : | deduites de ces critéres, un accord satisfaisant. 


| 
| 
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Fraction molawe en methano/ dans /e liquide _——— oo 
Fig. 13. Systéme methanol-dimethylformamide. Courbe eptenennoes 
"1 titre moléculaire d'un constituant dans le liquide 
log — = f (z,). = titre moléculaire d'un constituant dans la vapeur 

"2 = pression d'expérience 
= tension de vapeur d'un constituant & une tem- 
, Page aes pérature déterminée 
Les systémes binaires Benzéne-Diméthylforma- y = coefficient d’activité 
mide et Méthanol-Diméthylformamide ne sont | A et B = constantes des équations de Van Laar 


CONCLUSIONS 
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R. C. L. Boswortn : Heat Transfer Phenomena. Associ- 
ated General Publications, Pty. Ltd., Sydney. John Wiley 
& Sons Inc., New York. Pergamon Press Ltd., London, 
1952. £3. 10s. (Australian) or $6.00 or 42s. 


THs monograph forms an extension and addition to the 
opening paper read at a conference on Heat Transfer at 
Sydney in 1948. As the author says in his preface, the 
book comprises “ essentially a review of recent develop- 
ments in Heat Transfer and is concerned with : the con- 
duction of heat in gases, radiative transfer, conduction in 
solids and liquids, forced and natural convection, heat 
flow with simultaneous phase change and the use of elec- 
trical and other model methods.” Indeed, the book is 
not meant as a handbook for those concerned with thermo- 
technical calculations, but merely provides an exposition 
of the present fundamental insight into the various fields 
of heat transfer as well as of the difficulties still to be 
overcome. In this respect the book is a very welcome 
addition to the well-known handbooks and enables the 
reader, with little effort, to become familiar with this 
more fundamental background of heat transfer problems. 


Owing to the size of the book, every problem is dealt 
with rather summarily. However, thanks to the lucidity 
of the author, who consistently treats all forms of heat 
transfer from the point of view of “ carrier-mean free 
path,”’ the work has become a most attractive exposition 
which reads like a novel for those who are at all acquainted 
with heat transfer problems. The numerous references 
at the end of every chapter present sufficient points of 
departure for a more profound study of the various 
subjects. 

However, the book shows that the author — who has 
already made several contributions to the literature in 
this field — evinces a partiality for some of the problems 
treated. This accounts for the fact that the book suffers 
from a certain lack of balance as regards the attention 
given to the various subjects. 


For example, in addition to the brief discussion of 
dimensional analysis in the chapters on natural and 
forced convection, 15 pages of the last chapter are devoted 
to a rather fruitless treatment of thermodynamical 
similarity, whereas the important subject of analogy be- 
tween heat and mass transfer is disposed of in slightly 
over 2 pages. “The use of the equivalent electrical 
circuit "’ receives 22 pages, but no attention is given to 
the recent development in the field of heat transfer at 
very high gas velocities. 

Apart from this criticism in comparatively minor points, 
the author may be said to have succeeded in treating the 
comprehensive subject of heat transfer in a concise but 
nevertheless extremely lucid way. In consequence the 
book is particularly suited for teaching purposes, although 
its rather high price prevents its extensive use as a text 
book. 





Finally, it would be commendable in a subsequent 
edition to expend as much care on the figures as has been 
bestowed on the printed text, and to consistently use the 
kg.cm.sec. Joule system instead of expressing the various 
quantities in arbitrarily selected units (radiative constant 
oe in erg.cm-?sec1k* and B.T.U.’s ft-? hr? °R4, 
radiancy of gases E in kgcal m~ hr-!, thermal conductivity 
k in cal. em-!, sec™!, °C], heat transfer coefficient A in 
B.T.U.’s ft? hr-! °F-!). This alteration would be wel- 
comed, especially in this book which forms a link between 
the technical and the fundamental science. 


C. v. HEERDEN. 


E. A. GuGGennem : Mixtures, Oxford University Press 


1952, 42s. 

L’auteur dont le traité * Statistical Thermodynamics,” 
rédigé en collaboration avec Fow er, est un classique 
dans tous les pays et dont le plus récent Thermodynamics 
est un exemple de précision, présente avec clarté et 
simplicité sa théorie des solutions strictement réguliéres. 
Aprés un rappel de thermodynamique classique et une 
introduction aux méthodes de la mécanique statistique, 
il expose la théorie des solutions liquides “ s. réguliéres ” 
(premiére approximation et approximations supérieures) 
qu'il met en paralléle avec la méthode quasichimique de 
Berne et celle de Kirnkwoop. I! compare ensuite les 
résultats théoriques et expérimentaux. 


De l'étude des solutions liquides diluées, lauteur 
passe a celle des imperfections dans les cristaux. Il 
étudie ensuite les réseaux cristallins surimposés (point A). 
Les mémes approximations que précédemment sont 
utilisées. L’auteur passe ensuite a la mécanique statistique 
des mélanges gazeux parfaits et réels en liaison avec la 
théorie des “ états correspondants.’”’ Un chapitre est 
consacré & la thermodynamique et a la statistique des 
solutions superficielles, nouvelle application de la théorie 
des solutions “s. réguliéres.”’ 


L’étude des solutions athermes et non athermes 
(molécules de dimensions différentes) et des solutions de 
macromolécules compléte cet exposé des théories simples 
particuliérement fécondes et dont les résultats peuvent 
souvent étre raisonnablement confrontés avec les 
expériences. 

Bien entendu, dés l'introduction, l’auteur insiste 
sur le fait que ces modeéles ne peuvent prétendre représenter 
que les solutions les plus simples. Mais ceci permet 
néanmoins de traiter quantitativement des problémes 
qui sont d’ailleurs suffisamment compliqués pour étre 
intéressants. 

A. Porier 
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J. H. pe Boer: The Dynamical Character of 
Absorption. 1953, O.U.P., 30/-, 239 p., XV. p.p., 45 figs. 
To his well-known work on the emission and absorption 
of electrons Dr. pe Borr has added this lively book 
on the physical adsorption of molecules. The importance 
of a mechanical picture of the kinetic processes of physical 
adsorption and of the conditions in layers of molecules 
is strongly emphasized. Because they have been the 
subject of an article by the author (in Advances in Catalysis 
Volume II) energy relationships do not receive prominence. 


The introduction of the book outlines the kinetic 
approach: the terms in the fundamental equation 


Collision frequency _ Lifetime on Number of 
x = 
with surface surface adsorbed molecules 


are dealt with individually in the first few chapters 
and it is here that the pictorial aspect is stressed. Perhaps 
overstressed ; for a discussion of energy and entropy of 
adsorption we have to be content with brief accounts 
in later chapters. 


The latter and larger part of the book considers the 
properties of the adsorbed layer in terms of a 2-dimensional 
gas. Beginning with ideal behaviour this account goes 
on to introduce constants (b, and a,) corresponding to 
those of the van der Waals equation for the 3-dimensional 
gas. These constants form the basis of a treatment of 
condensation :nd critical phenomena in two dimensions 
which is extended to adsorption equilibrium — monolayer, 
multilayer and finally condensation. The last chapter 
of the book is virtually an appendix on capillary effects. 


The opening chapters have an attractive vigour and 
clarity which will appeal to most readers. The clarity 
of later sections might have been increased further if 
chapters had been here added on the energetics of adsorp- 
tion and the thermodynamic prcperties of mobile adsorbed 
layers. The use of van de Waals equation which has tne 
advantage of familiarity could have been supplemented 
by a clear summary of relations between it and its 
constants with those of a virial expansion. 


Despite this unevenness the book is a unique and 
stimulating account of physical adsorption, and will 
be of great value to all who deal with these problems. 

P. Gray. 


F. Witp : Estimation of Organic Compounds, 239 p.p. 
Cambridge University Press, 1953, 25s. net. 
The author describes different methods for the estimation 
of groups in organic compounds and of certain compounds, 
e.g. the olefines, hydroxyl groups containing compounds, 
mercaptans, carbonyl compounds, amines, nitro, cyano 
and nitroso compounds. 

The reviewer thinks that Dr. Wiip has been succesful 
in the choice of the compounds in question, many of these 





being of importance for the organic chemical industry. 
The description of the methods is thorough. 


In determining the hydrogen value of olefines it is 
useful to investigate the hydrogenated product, in order 
to check if it is really the expected product. It is known 
that for the examination of oils and fats the thiocyanogen 
number is a “ distinctive number” having a certain 
significance, only when all the conditions during the 
determination are fixed within very definite limits. 


Stress could be laid still more on the necessity for 
the application of the hydrogen values. The author does 
this for gaseous olefines ; he says : “* catalytic hydrogena- 
tion is the only rapid method, giving a reliable indication 
of the degree of unsaturation.” 


The halogen numbers are often not correct, the halogen 
addition generally depends on an equilibrium. 


VAN Der Streur showed that the equilibrium constants 
are dependant not only on the temperature (they rise with 
decreasing temperature), but also on the nature and the 
position of the double bond. 


This book is very valuable, it can be recommended 
to the research and analytical laboratories, which have 
to investigate or to synthetize the compounds discussed 
in the book. It has an excellent subject and author index. 

H. I, Warerman 


Food Investigation, 1952. Published by Dept. of 
Scient. and Industr. Res. H.M. Stationery Office, 1953. 
56 pp. Price 2s. (U.S.A.: 50 cents). 

This report covers the wide range of activities of the Food 
Investigation Board, and gives an interesting general 
picture of the problems which are being tackled. A large 
proportion of the work is concerned with the development 
of improved methods of storage, particularly of meat, fish 
and fruit and vegetables, so as to prevent deterioration of 
quality. An interesting point is that, whereas some fish 
is found to have deteriorated after a month in cold store, 
its quality improves again after a further two months’ 
storage. 

Chemical Engineers will be interested in the work on 
heat transfer at low temperatures. Research is being 
carried out in an air blast tunnel freezer, on the effect of 
air temperature, velocity, turbulence and relative humidity 
on the rate of heat transfer to slabs of fish ; temperatures 
are measured by means of thermocouples in the solid and 
in the air stream. Freeze concentration of beer and 
evaporation at low temperatures are being studied in 
co-operation with the brewers’ organizations. Improved 
methods of drying fruit and vegetables and the dehydra- 
ting of meat are being developed. An incidental problem 
is the tin plating of aluminium for use in canning, in order 
to prevent corrosion and to enable it to be soldered. 


J. F. Ricwarrson, 
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A. G. Gayvon and H. G. Wotrnarp : Flames: Their 
structure and radiation. Chapman and Hall, 1953. 
340 pages. Price 55s. 


A great deal of work has been done on flames since Bone 
and Townsenp’s Flame and Combustion in Gases was 
published in 1927, and numerous texts have appeared 
since then. Outstanding among those available in English 
today are the translation of Jost’s excellent Explosion 
and Combustion Processes in Gases (1935) and Lewis and 
von Ex.se’s Combustion Flames and Explosions in Gases 
(1951). These two texts have become standard works for 
students of flame. To them has now been added Flames 
by Dr. Gaypow and Dr. Wotruarp. The authors’ close 
connection with nearly every aspect of the subject they 
discuss is obvious. It is a valuable contribution. 


The early chapters (I, I1) are concerned with flames in 
premixed gases. A general descriptive account is given of 
their properties ; more general material may also be found 
in a later chapter (XIV). A detailed account (111) of flow 
patterns in flames and of the shapes of flame fronts 
follows : the influence of sound on flame shape and propa- 
gation, and the accoustic effects of flames are not con- 
sidered until a later chapter (VII). Experimental methods 
of measuring flame speed are listed and their reliability 
discussed ; this section (IV) receives particular emphasis. 
The mechanism of flame propagation is then considered 
(V) and a discussion of turbulent flames is included ; 
chemical] kinetics of combustion processes are not discussed 
though the mechanism of carbon formation is considered 
in some detail (VIII). Both thermal and diffusion theories 
of steady flame propagation are considered and various 
explosions are quoted relating burning velocity to the 
other parameters typical of the flame and the chemical 
reactions occurring in it. The authors consider tnat 
problems of spontaneous ignition, critical ignition energy 
and flame stability, which might have found a place here, 
and which have been dealt with at length by Jost, and 
by Lewis and von Ese, to be beyond the scope of this 
book. They also deliberately avoid purely mathematical 
discussion ; but a more critical comparison and assessment 
of the basic assumptions of different theories would have 
been a valuable contribution to this section. Diffusion 
flames are described in chapter VI. The results obtained 
by the authors by the spectroscopic examination of flat 
diffusion flames are described, and diffusion and pre-mixed 
flames are contrasted. Turbulent diffusion flames, furnace 
flames, flames of droplets, of fuels on wicks and of the 
candle are considered briefly. Solid particle combustion 
is mentioned later (XIV). A chapter (IX) follows on 
radiation processes in flames. This section centres on the 
problems of dis-equilibrium both of energy distribution 
and of chemical reaction and forms a valuable intro- 
duction to the discussion of flame temperature (X, XI, 
XII). The validity of the concept of temperature, and 
the theoretical and practical limitations of the various 
methods of measurement are discussed with clarity. This 





valuable section ends with the systematic calculation of 
flame temperatures ; typical results are given for a variety 
of flames. Ionization in flames (X U1), a process notorious- 
ly out of equilibrium, follows aptly the radiation and 
temperature sections and ends the formal topic of the 
book. Chapter XIV “Some unusual flames” which 
follows, contains a miscellany of recent information ; 
various unresolved problems are summarized in Chapter 
XV. 


Some lack of logical development may be felt in the 
authors’ exposition. But no one exclusive arrangement 
of subject matter exists and its lively treatment here 
affords ample compensation. It is an excellent book. 


Perer Gray. 


Cart Mirrac: Die Hartzerkleinerung (The Milling of 
Hard Materials) Springer-Verlag, Berlin/Géttingen/ 
Heidelberg, 1953. vi + 342 pages, 190 illustrations. 
DM 40.50. 


Cart Mirrac is the Grand Old Man of the practice of 
milling hard materials in Germany; for more than 50 
years he has worked in this field, which has expanded so 
rapidly in the last decades that a comprehensive survey 
would be a monumental task. Murrac therefore selects 
his material be to of use to students in technical colleges 
and mining schools, to constructional and chemical 
engineers, and to industrial technologists and scientists, 
who may wish to learn the principles of milling hard 
materials, to choose the right mill for specific jobs or to 
improve upon existing types of mill. 


The book is arranged in three parts. The first, after a 
short historical introduction, describes 27 types of crushing 
and milling machines. The second deals with the funda- 
mental physics of the milling operation and the third 
describes ancillary equipment and experience in milling 
a number of important materials. 


The first part, amounting to half the book, is most 
attractively arranged in a consistent pattern for each 
type of mill: an introduction, a description of the plant 


and its field of application. Of the 190 illustrations, all 
apparently drawn specially for the book, 129 are in the 
first part. Most of them are sections and show very well 
how the mills work, an advantage to readers with little 
or no German ; however, only the outside of the Hardinge 
mill is drawn, possibly because the author says on p. 150 
that the principle of its working is unsound ; the reader 
will want to know .why. Apart from a few apparent 
prejudices, this part of the book should indeed be useful 
to the classes of people to whom it is offered. By the use 
of numerous formulae, tables and graphs, it is possible, 
for many of the more important minerals of commerce, 
to determine the size of mill for a given output and grist. 





Book reviews 


The theoretical part deals with some aspects of the 
fundamentals of comminution and gives some practical 
hints Mrrrac has found useful. Chiefly following the 
work of von Rirrincen and Kick, and of Smexkar and 
ANSLEM, Mrirrac gives a useful summary of the relation- 
ship between the work done to produce new surface area 
in milling ; the * efficiency’ of a mill can be judged by 
his calculations of this relationship, expressed throughout 


2 
as 50 many gcm/cm”. 


For the measurement of particle size only ANDREASEN’s 
pipette is mentioned, and the pages on micromeritics are 
curiously old-fashioned. The reason for this can be seen 
by referring to the list of references at the end of the 
book ; 21 are pre-1928, 55 are 1928-37, and only 14 are 
post-1937 (up to 1949). None of the 19 English and U.S. 
papers referred to is after 1936, and 6 out of the 19 have 
spelling errors, including Srock’s law (also in the text) 
and A. F. TaGccerr. 


Both in parts one and two a lot of information is given 
on the right speed and load of ball mills, but CoGHti. 
and pe VaANry’'s and others’ useful formulae for calculating 
the right size of balls is not given, much less OLEvsky’'s 
more recent formula (1948). 


Part two is concluded with a section on future develop- 
ments in milling plant techniques based on explo_ion, 


the * puffed wheat’ and vibration principles, but some- 





thing might have been said of ‘ micronising,’ fluid energy 
mills and ultrasonic comminution. 


The third part is intended to help the reader to find 
‘the right machine for the right job,’ and has a valuable 
section on ancillary plant such as feeding and classifying 
plant, including sieving machines, wind-sifters, and dust 
separators. Wet classifying is illustrated by only one type 
of plant, but the other types of plant are as well illustrated 
and described as in part one. 


The choice of plant for a given raw material is illustrated 
by descriptions of typical processes for the comminution 
of iron, copper, lead-zinc, tungsten, and tin ores ; of coal, 
though not a hard material ; of carbonate rocks, feldspar 
and quartzite ; of cement slags, etc. ; of calcium carbide 
and cyanamide ; of phosphates and of carbon for making 
electrodes. 


Car Mirtac does not deal very much with the problems 
of extremely fine grinding, possibly because this is usually 


considered as a separate technical category in Germany 
( Feinstzermahlung) ; but few men can have had so much 
experience in the classical methods of milling in the heavy 
chemical industry, and his descriptions of the different 
sorts of mill will continue to be useful for many years 
because the older types, yes even millstones, will be with 
us for centuries. 

R. H. 5. Ropertson. 





